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SUMMARY 


The  objective  of  Project  THEMIS,  A  Center  for  Remote  Sensing 
(Contract  No.  DAAK02-68-C-0089)  at  the  University  of  Kansas  is  advance¬ 
ment  in  the  fields  of  remote  sensing  with  primary  emphasis  on  user  appli¬ 
cation  of  remotely  sensed  data  of  natural  and  cultural  environments. 

The  program  includes  four  major  areas  of  interest:  theoretical  studies,  system 
design,  data  processing  and  display,  and  the  data  analysis  and  interpreta¬ 
tion.  During  this  report  period  work  has  been  performed  along  the  lines 
drawn  in  the  study  plan  and  is  a  continuation  of  the  work  reported  in  the 
ureyious  semi-annual  technical  report  under  this  contract.  Theoretical 
investigations  of  the  frequency  dependence  of  ultrasonic  scatter  from 
statistically  known  rough  surfaces  and  the  results  of  the  ultrasonic  model 
studies  experiment  are  reported.  The  design  and  construction  of  the  one 
octave  (4-8  GHz)  broadband  radar  system  is  almost  complete;  it  is  presently 
in  the  final  testing  and  calibration  phase.  Except  for  a  few  delays  in  de¬ 
livery  of  equipment  placed  on  order  for  the  interface  units  between  the  exist¬ 
ing  analog  IDECS  and  digital  equipment,  the  progress  in  the  construction  of 
the  data  processing  and  display  system  is  as  anticipated  in  the  study  plan. 

As  described  in  the  data  processing  section  of  this  report,  studies  are  reported 
on  the  application  ox  adaptive  classification  techniques  to  both  spatial  and 
measurement  space  clustering  and  on  the  application  of  Bayes*  statistical 
decision  theory  to  data  obtained  from  a  remotely  sensed  environment. 

Project  THEMIS  at  the  University  of  Kansas  is  funded  under  ARPA-DOD 
Work  Order  1079  and  is  monitored  as  Contract  DAAK02-68-C-0089  by  the 
U.S.  Army  Engineer  Topographic  Laboratories,  the  Geographic  Sciences 
Division.  Mr.  Bernard  B.  Scheps  has  been  the  Project  Engineer  since 
contract  inception;  Mr.  Abraham  Anson  has  recently  replaced  Mr.  Scheps 
in  this  capacity.  The  project  is  supported  at  the  University  of  Kansas  by 
the  following  personnel,  some  of  whom  are  not  directly  funded  under  THEMIS, 
but  whose  work  is  contributory  to  the  effort: 


A.  Faculty  Members 

Electrical  Engineering  6 

Geography  1 

Geology  2 


B.  Ph.D.  Candidates 

Electrical  Engineering  15 

Geography  5 

Geology  2 


C.  Masters  Candidates 

Electrical  Engineering  7 

Geography  5 

Geology  2 


In  addition  to  the  personnel  listed,  a  number  of  undergraduate  students 
are  employed  as  technicians  under  the  three  principal  disciplines.  A 
working  arrangement  is  being  developed  with  faculty  members  of  the 
departments  of  Psychology  and  Sociology  in  the  field  of  pattern  recognition 
and  false-color  perception  and  some  experimentation  was  carried  on  during 
this  report  period.  It  is  hoped  that  a  formalized  arrangement  can  be  made  in 
the  future  to  carry  this  work  forward  in  a  systematic  manner. 

Since  its  inception  in  1964,  the  Remote  Sensing  Laboratory  at  the 

University  of  Kansas  has  been  actively  engaged  in  a  number  of  studies 

which  relate  directly  to  the  work  being  done  under  Project  THEMIS.  Themis 

is  of  particular  importance  to  the  University  of  Kansas  since,  because  of 

its  broad  objectives,  it  provides  a  focus  for  the  work  being  done  under  the 

other  studies.  As  a  result,  considerable  cross -fertilization  between  pro- 

* 

Jects  is  occuring  almost  automatically  and  the  total  effort  of  the  Remote 
Sensing  Laboratory  is  emerging  as  a  cohesive  entity. 

The  studies  currently  underway  at  the  Remote  Sensing  Laboratory 
in  addition  to  Project  THEMIS  are: 


Name  of  Project 


Granting  Agency 


No. 


Theoretical  and  Experimental 
Study  of  the  Nature  of  Rough - 
Surface  Scattering 

♦Research  in  Radar  Scattero- 
metry  and  Altimetry  Applied 
to  Oceanographic  Problems 

Study  of  Waves  Backscattered 
from  Layered  Rough  Surfaces 

♦Study  the  Utility  of  Radar  and 
Other  Remote  Sensor  Imagery  in 
Thematic  Land  Use  Mapping  for 
Geographic  Research 

♦Radar  Studies  Related  to 
Orbiting  Research  Labs 

A  Study  of  High  Performance 
Ai  ienna  Systems  for  Deep 
Sp  ce  Communication 

♦Multi-Image  Correlation 
Systems  Study  for  Military 
Geographic  Intelligence 

Application  of  Nonlinear  Opti¬ 
mization  Techniques  to  Obtain 
Network  Designs  having  Low 
Multiparameter  Sensitivity 

Study  of  Backscattering 
Characteristics  of  Mars 
and  Venus 


NSF 

USNOO 

NSF 

USGS 

N  ASA/MS  C 
NASA/MSC 

USAETL 

NSF 


Ryan  Aero¬ 
nautical  Co. 


A  Study  of  Optimum  Interpre-  NASA 

tability  of  Digitized  Image’y 

♦Project  Darien  USAETL 

♦Project  SAND  USAETL 


★ 


Geoscience  Participation  or  Geoscience  Project 


GK  1153 

NG2306-67C-0044 

GK  875 

14-08-0001-10848 

NAS  9-7175 

NGR  17-004-013 

DAAK02-67-C-0435 

GK  3657 

BG1 7-002-053 

DAAK02-68-M-6993 

DAAK02-69-M-0255 


PROPOSALS 


♦Research  Contract  to  U.S.  Dept.  ofAgric. 

Study  the  Utility  of  Radar  in 

Applications  and  Research  on 

Crops,  Natural  Vegetation, 

and  Soil  Moisture 

*A  Center  for  Excellence  in  NSF 

Remote  Sensing 

Theoretical  and  Experimental  NSF 

Study  of  the  Nature  of  Rough 
Scattering — A  Continuation 


PROJECT  THEMIS:  A  CENTER  FOR  REMOTE  SENSING 
SEMI-ANNUAL  TECHNICAL  REPORT  (U) 


I.  INTRODUCTION 

The  purpose  of  this  report  is  to  describe  the  technical  work 
accomplished  under  contract  No.  DAAK02-68-C-0089  during  the  period 
April  1,  1968  through  October  1,  1968.  It  is,  in  effect,  a  continuation 
of  the  work  described  in  the  previous  semi-annual  technical  report^ 
with  emphasis  on  the  rate  at  which  the  work  under  the  various  tasks  has 
been  progressing  and  the  technical  and/or  scientific  significance  of  the 
results  obtained  to  date.  This  report  has  been  prepared  in  the  form  of 
a  summary  in  which  appropriate  references  are  made  to  the  technical 
reports  and  memoranda  which  describe  the  work  under  each  subtask  in 
detail.  Some  of  the  more  significant  memoranda  and  parts  of  long  techni¬ 
cal  reports  are  appended  in  this  text. 

Except  for  a  few  delays  in  delivery  of  equipment  placed  on  order, 
the  range  of  progress  of  the  overall  project  is  as  anticipated. 


II.  TECHNICAL  REPORT 

TASK  4.1*  TARGET -SENSOR  INTERACTION  STUDIES 

This  task  includes  theoretical,  modeling,  and  experimental  studies 
of  broad-spectrum  radar  systems.  Considerable  progress  has  seen  made  in 
all  of  the  above  studies,  as  detailed  below. 


*Task  numbers  are  those  of  the  USAETL  Purchase  Description  dated  16  September 
1967. 


SUBTASK  4.1.1  THEORETICAL  STUDIES 


Theoretical  studies  of  broad -spectrum  scattering  are  being  performed 
at  the  University  of  Kansas  under  National  Science  Foundation  grants. 

Major  work  under  the  above  heading  appeared  as  a  technical  publication 
in  the  1968  Wescon  convention  of  the  IEEE  .  This  publication  is  appended 
as  Appendix  A. 

SUBTASK  4.1.2  ULTRASONIC  MODEL  STUDIES 

Technical  Memorandum  113-2  (March  1968),  which  appeared  under 
Appendix  A  of  the  previous  semi-annual  technical  report  of  Project  THEMIS ^  , 
examined  the  v?  ious  aspects  involved  in  the  study  of  the  frequency  dependence 
of  backscatter  from  statistically  rough  surfaces  using  a  broad  spectrum 
backscatter  acoustic  system.  This  included  a  brief  clscussion  of  the 
acoustic  simulation  technique,  a  detailed  description  of  the  system  instru¬ 
mentation  and  calibration  problems,  and  the  data  processing  and  analysis 
involved  in  the  scattering  and  image  experiments.  Since  March  1968,  the 
program  has  progressed  through  the  final  stages  of  data  acquisition  and  analysis. 

The  complete  work  is  compiled  and  presented  by  Dr.  J.  W.  Rouse,  Jr.  as 

(3) 

his  Ph.D.  dissertation'  ' .  An  interim  technical  report  under  Project  THEMIS 
in  the  form  of  Dr.  Rouse's  dissertation  is  presently  in  the  printing  process. 
Appendix  B  is  a  condensed  form  of  the  results  reported  in  the  above  disserta¬ 
tion. 

SUBTASKS  4 . 1 . 3 ,  4.1.4,  4.1.5,  4.1.6  ELECTROMAGNETIC  INSTRUMENTATION 
AND  MEASUREMENT 

These  subtasks  are  combined  into  a  single  task  description  because 
of  the  close  inter-relation  of  subtasks  which  prevents  convenient  separation. 


*NSF  Grants  GK-875  and  GK-1153 
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The  one  octave  bandwidth  system  being  constructed  is  near 
laboratory  operation.  This  system  was  operated  in  the  laboratory  in  August 
1968,  however,  incorrect  manufacturers  rating  led  to  the  burn-out  of  the 
PIN  diode  switch  used  for  duplexing.  This  malfunction  required  redesign 
of  the  switching  system  to  utilize  higher  power  rating  switches.  Hewlett 
Packard  Associates  is  furnishing  the  additional  replacement  switches  required 
at  no  cost  due  to  the  error  in  specification  which  led  to  the  initial  failure. 

The  replacement  switches  are  due  by  November  1  and  calibration 
and  alighment  of  the  RF  portion  of  the  system  can  be  started  with  the  new 
switching  arrangement. 

CRES  Technical  Memorandum  133-1-2,  shows  the  latest  detailed 
block  diagram  and  component  diagram  of  the  system  after  the  above  mentioned 
change  in  duplexing.  The  GFE  boom  truck  has  been  repaired  and  Is  now 
ready  for  limited  range  use.  The  original  bucket  has  been  replaced  by  a 
mounting  fixture  for  the  antenna  positioner.  The  antenna  positioner  itself 
is  the  only  remaining  item  not  yet  received;  however,  delivery  is  estimated 
by  November  1 .  The  delay  in  ordering  the  positioner  was  due  to  the  time 
required  for  fund  transfer  to  purchase  the  item  and  not  to  lack  of  definition. 

The  Logic  Control  system  required  to  perform  the  gating  of  the 
pulses  for  subsequent  integration  in  time  or  frequency  is  described  in 
CRES  Technical  Memorandum  133-1-3.  This  system  is  now  in  the  final 
construction  stages.  All  circuits  have  been  through  initial  breadboard  and 
test  and  the  svstem  is  being  fabricated  in  replaceable  plug-in  card  modules. 

The  GFE  laboratory  van  and  truck  have  been  received;  however,  this 
vehicle  has  not  yet  been  placed  in  operating  condition. 

CRES  Technical  Memoranda  133-1-2  and  133-1-3  are  appended  as 
Appendices  C  and  D. 


TASK  4.2  SENSOR  AND  PREPROCESSOR  DEVELOPMENT 

The  emphasis  in  this  task  is  on  parameter  and  design  considerations 
and  development  of  prototype  sensors  and  sensor  systems  including  associated 
preprocessing  of  sensor  data.  As  was  mentioned  in  the  previous  semi-annual 
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technical  report ,  due  to  fund  limitations  subtasks  4.2.2  (Cartographic  radar 
systems)  and  4.2.5  (semi-focused  imaging  systems)  are  not  being  attacked 
urder  this  contract.  Work  on  subtask  4.2.1  (evaluation  of  broad-spectrum 
imaging  systems)  is  still  pending  until  data  from  the  system  under  develop¬ 
ment  are  available . 

SUBTASK  4.2.3  STUDY  OF  IMAGING  RADAR  SYSTEMS  FOR  SMALL  SATELLITES 

At  this  time  work  is  being  done  on  defining  a  class  of  ground -mapping 
imaging  radar  systems  to  be  flown  aboard  small  earth  orbiting  satellites. 

This  study  is  being  performed  under  NASA  Contract  No.  NAS  9-7175.  Investi¬ 
gations  made  in  this  study  which  are  applicable  under  the  above  subtask 
title  will  be  reported  both  to  the  granting  agency  and  to  USAETL. 

SUBTASK  4.2.4  A  STUDY  OF  RADAR  SCATTEROMETER  SYSTEMS 

Study  of  radar  scatterometer  systems  is  presently  being  implemented 
Jointly  under  the  auspices  of  the  Naval  Ordnance  Laboratories  (Contract  No. 
N62306-67-C-0044) ,  NASA  (Contract  No.  NAS  9-7175)  and  Project  THEMIS. 
Parameter  studies  are  being  conducted  for  sea-ice  measurements  and  include 
considerations  of  roughness,  frequency,  resolution  area,  polarization  and 
antenna  patterns . 

(4) 

A  final  memorandum'  1  on  the  correlated  split-gate  altimeter  has  been 
completed  under  USNOO  Contract  NG2306-67C-0044 .  With  the  results  of  the 
range  error  statistics  for  square  law^  and  linear ^  detection  of  independent 
pulses,  the  correlated  or  dependent  pulse  analysis  provides  a  theoretical 
base  for  experimental  studies  on  the  split-gate  altimeter.  The  split -gate 
tracking  system  was  studied  because  it  is  frequently  employed  in  the  measure¬ 
ment  and  tracking  of  altitude  with  pulsed  radar  altimeters .  The  accuracy  is 
affected  by  the  deviation  caused  by  random  fluctuation  of  returns  from  surface 
scatterers .  The  resultant  jitter  creates  range  or  altitude  noise  in  altimeter 
response . 

(71 

Analysis'  '  of  the  backscattering  cross  section  for  sea  states 
generated  by  fully  developed  winds  has  been  completed  in  October,  1968. 
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The  empirical  sea  wave  spectra  for  different  winds  has  been  combined  with 
rough  surface  scattering  theory  to  provide  angular  dependence  of  cross 
section  (aQ  versus  theta)  for  these  winds. 

A  series  of  scatterometer  measurements  made  in  the  Arctic  Ocean 
north  of  Point  Barrow,  Alaska,  have  demonstrated  the  capability  of  the 
13.3  GHz  system  to  differentiate  different  sea  ice  categories^}  Differ¬ 
entiations  have  been  achieved  in  winter  and  polar  ice  types . 

TASK  4.3  DATA  PROCESSING  AND  DISPLAY 

The  purpose  of  this  task  is  to  improve  existing  techniques  for  data 
display,  such  as  those  employed  in  the  IDECS  system,  and  to  originate  new 
techniques  for  processing  and  display  of  remote  sensor  data  for  improved 
interpretability.  At  this  time,  the  working  subtasks  under  4.3  fall  into 
two  categories,  1)  improvement  of  the  IDECS  system  (Subtasks  4.3.1  and 
4.3.3)  and  2)  implementation  of  adaptive  and  Bayes  decision  techniques 
(Subtasks  4.3.4  and  4.3.5).  As  work  progresses  in  this  categories, 
subtasks  4.3.2  (evaluation  of  color  combinations  of  multiple  images)  and 
4.3.6  (new  display  techniques)  will  be  implemented.  Ultimately,  the 
adaptive  and  Bayes  decision  techniques  will  be  instrumented  and  inter¬ 
faced  with  the  IDECS  system  in  such  a  way  that  the  IDECS  operator  will 
have  maximum  flexibility  in  enhancing  data  by  direct  combinations  of  images 
along  with  computer  processed  data. 

Representative  of  the  work  accomplished  to  date  under  task  4.3 
are  thiec  memoranda  which  are  appended  as  Appendices  E,  F  and  G.  The 
first  of  these  (Appendix  F)  deals  with  the  application  of  adaptive  classifica¬ 
tion  techniques  to  both  spatial  and  measurement  space  clustering.  The 
application  of  Bayes'  statistical  decision  theory  to  data  obtained  from  a 
remotely  sensed  environment  is  considered  in  Appendix  F.  The  progress 
and  mechanization  of  the  IDECS  system  improvements  are  outlined  in  Appendix 
G. 


TASK  4.4  DATA  ANALYSIS  AND  APPLICATION 


This  task  involves  in  part  the  analysis  of  the  data  obtained  from  the 
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broad -spectrum  radar  systems  of  task  4.1 .  Also  included  are  studies  relating 
to  the  increase  of  information  content  in  multiple  radar  images  over  single 
Image  content. 

Subtasks  4.4.2  (relation  of  image  data  to  wavelength,  polarization, 
depression  angle,  resolution,  and  bandwidth)  and  4.4.3  (information  increase 
with  the  use  of  multiple  radar  images)  will  not  be  implemented  until  broad- 
spectrum  image  data  are  available. 

SUBTASK  4.4.1  EVALUATE  GEOGRAPHIC  AND  GEOLOGIC  INFORMATION 
AVAILABLE  FROM  SINGLE  RADAR  IMAGES 

Since  no  data  will  be  available  from  the  broad-spectrum  radar 
systems  of  task  4.1  before  the  second  year,  the  efforts  under  this  task 
are  devoted  to  geographic  and  geologic  evaluations  of  single  image  radar 
data.  During  the  second  and  succeeding  years  this  task  will  be  changed 
to  an  evaluation  of  the  IDECS  and  associated  signature  selection  system 
as  aids  in  image  interpretation.  A  final  report  by  Drs.  Peterson  and  Simonett 
on  the  geomorphology  of  the  Wasatch  Range  (Appended  as  a  preliminary 
report.  Appendix  F,  in  the  previous  report)  will  be  issued  as  a  special 
Interim  report  in  the  near  future . 

Work  on  single  radar  images  has  been  done  at  CRES  under  Project 

SAND  (Contract  No.  DAAK02-69-M-0255)  to  study  the  availability  of  construc¬ 

ts) 

tion  materials  in  the  Mekong  Delta'  ' .  If  more  information  on  this  work  is 
desired,  the  reader  can  be  supplied  with  a  copy  of  CRES  memorandum  156-1 
which  describes  the  investigations  under  Project  SAND  to  date. 


III.  CONCLUSIONS  AND  RECOMMENDATIONS 

Except  for  a  few  delays  in  delivery  of  components  placed  on  order  for 
subtasks  4.1.3  -  4.1.6  and  task  4.3,  the  rate  of  progress  of  the  overall 
program  of  Project  THEMIS  is  as  anticipated  in  the  study  plan.  Due  to 
these  delays,  subtasks  4.1.3  -  4.1.6  are  about  90  days  behind  schedule;  It 


*DAAK02-1560-824 


is  anticipated  that  the  broadband  radar  system  will  be  in  operation  by  early 
December  1968 . 

Appendices  A  through  G  of  this  report  present  the  reader  with  a  detailed 
picture  of  the  efforts  made  under  Project  THEMIS  since  April  1968.  Most  of 
the  work  was  concentrated  under  tasks  4.1  and  4.3.  The  output  under 
task  4.2  was  limited  by  shortage  of  funds,  and  major  work  on  task  4.4  will  be 
commenced  after  broad -spectrum  image  data  are  available  and  the  IDECS 
system  is  fully  operational . 
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SUMMARY 


FREQUENCY  DEPENDENCE 


An  acoustic  simulator  Is  used  to  model 
radar  backscatter  and  the  frequency  dependence 
of  the  return  signals  from  statistically  known 
rough  surfaces  Is  studied.  One  of  the  model 
statistically  rough  surfaces  Is  made  of  rubber 
polymer  obtained  with  a  clay  mold  and  the  other 
is  made  of  steel.  The  profiles  of  the  surfaces 
are  measured  with  p  depth  gauge  so  that  the  cor¬ 
relation  functions  and  the  standard  deviations  of 
the  surfaces  could  be  calculated. 

/* 

On  each  surface  average  backscattered 
power  are  measured  at  four  different  frequencies 
ranging  from  0.5  MHz  to  1.5  MHz.  Comparison 
with  the  scattering  jaw  obtained  under  Klrchhoff's 
approximation  s  hpws  real  good  agreement  when 
thp  parameters  in  the  scattering  law  are  properly 
Interpreted.  It  Is  believed  that  In  the  past  the 
scattering  law  could  not  be  verified  due  to  an 
Incorrect  interpretation  of  one  of  Its  parameters . 
Certain  precautions  that  must  be  taken  in  inter¬ 
pretating  the  scattering  law  are  also  pointed  out. 

INTRODUCTION 

This  study  concerns  the  change  in  angu¬ 
lar  variation  of  the  average  backscattered  power 
from  statistically  rough  surfaces  as  a  function  of 
the  Incident  frequency.  It  Is  clear  that  the 
angular  variation  of  the  return  signal  will  depend 
on  both  the  Incident  frequency  and  the  nature  of 
the  surface  under  study.  Fora  theoretical  pre¬ 
diction  assumptions  must,  therefore,  be  made 
about  either  the  correlation  function  of  the 
surface  or  its  Bessel  transform,  the  roughness 
spectrum.  A  case  of  great  practical  Interest  Is 
when  the  correlation  function  Is  monotone  de¬ 
creasing  but  may  have  oscillating  decaying 
tails.  For  simplicity,  we  assume  that  the  tail 
is  small  and  that  the  surface  height  distribution 
is  Gaussian.  We  shall  first  establish  some 
theoretical  predictions  for  two  different  cases  of 
surface  roughness  so  as  to  obtain  a  qualitative 
picture  of  the  backscatterlng  behavior  as  a  func¬ 
tion  of  frequency.  These  predictions  will  then 
be  compared  with  experimental  results.  Certain 
difficulties  in  interpretating  experimental  results 
are  pointed  out.  It  Is  found  that  previous  Inter¬ 
pretation  of  the  para  meters,  that  appear  In  the 
scattering  law  obtained  under  Klrchhoff's  approx¬ 
imation,  Is  Incorrect  and  a  new  Interpretation  is 
presented. 


For  a  plane  wave  Incident  upon  a  random 
surface  with  Gaussian  height  distribution, 
Klrchhoff's  method  predicts  an  average  far  zone 
backscattered  power  of  the  form* 

P(0)  *=  - — y-  \lo(2ksln0C)e‘K(1‘r(e))Cdc  (1) 

(  XcosQ)  J 

where  K  *  4k202  cos^O;  0  is  the  Incident 
angle;  k  ■  2r/X,  the  Incident  wave  number;  o» 
the  standard  deviation  of  the  surface,  r({)  the 
correlation  function  of  the  surface,  and  JQ  Is  the 
zero  order  Bessel  function.  In  (1)  we  have  left 
out  a  proportionality  constant  unimportant  for  our 
discussion.  The  integration  in  (1)  Is  over  the 
illuminated  area. 

If  X  is  large  compared  with  0 ,  so  that 
K  is  small  compared  with  unity,  (1)  may  be  ap¬ 
proximated  as  follows 

HO)  -  - - - j  e"K  \jn(2ksln00ede  + 

(Xcoser  J  0 

+(l/(XC0se)2]  Ke"K  lr(5)  Jo(2ksin00ed{«= 
■=P,(9)  +  P2(0)  J  (2) 

For  larqe  enough  Illuminated  area  the  first  term 
In  (2)  contributes  only  at  0=0.  The  scatter¬ 
ing  component  Is  the  second  torn  which  may  be 
rewritten  in  terms  of  the  roughness  spectrum  of 
the  surface,  W, 

P2(P)  ■=  (1  A)4  16*2  a2  e‘KW(2ksln0)  (3) 


Eq.  (3)  Is  In  agreement  with  the  first 
order  result  obtained  by  small  perturbation 
method.2  The  results,  however,  are  not  identi¬ 
cal  because  (3)  takes  Into  account  higher  order 
surface  slope  terms.  Looking  at  (3)  we  are  tempt¬ 
ed  to  say  that  It  predicts  an  angular  variation  of 
the  form  exp  (-K)  W  (2k  sin  0).  Actually,  things 
are  not  as  straight  forward,  since  we  always 
have  finite  beamwldth  In  performing  our  experi¬ 
ments.  Thus,  If  our  beamwldth  Is  0O  and  the 
incident  angle  is  10*,  we  cannot  use 
W  (2k  sin  10*)  but  rather  some  average  of 
W  ( 2k  sin  (10*+  9q/2)  ]  .  Hence,  at  a  fixed 
Incident  angle,  "we  are  not  concerned  with  Just  a 
single  k'  in  W(k'),  but  a  range  of  values  of  k'. 
Unless  we  maintain  the  same  beamwldth  as  we 
change  frequency,  the  problem  will  be  further 
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complicated,  since  2k  Increases  with  frequency 
and  our  beamwldth  Is  likely  to  decrease  with 
frequency.  Similar  kind  of  care  must  also  be 
exercised  when  we  consider  frequency  effect  on 
return  power  at  a  fixed  angle.  Thus,  It  Is  clear 
that  even  If  we  know  W(k*)<  a  precise  theoreti¬ 
cal  prediction  Is  still  difficult.  One  general 
conclusion  for  the  case,  X  »o.  Is  that  as  the 
Incident  frequency  Increases  the  angular  varia¬ 
tion  In  P2  (0)  becomes  smoother  when  the  beam- 
width  Is  kept  constant.  For  a  monotone  decreas¬ 
ing  W(k')  we  expect,  therefore,  P2  (0 )  to  be¬ 
come  flatter  as  frequency  Increases.  The  rate  of 
flattening  Is  smaller  for  a  W(k')  with  more  high 
frequency  components. 

Now  let  us  turn  to  the  more  practical  case 
when  K  Is  large  compared  with  unity.  Just  how 
large  K  should  be  has  never  been  established. 
Examining  (1)  shows  that  we  can  obtain  a  good 
approximation  by  applying  a  modified  saddle 
point  method.  In  the  past1,3,  the  practice  Is  to 
expand  r(£)  about  5=0,  and  take  the  first 
two  terms.  There  are  two  objections  to  doing 
this.  First  of  all,  5  =  0  Is  really  not  a  saddle 
point,  since  the  Integrand  vanishes  at  5=0. 
Seqondly,  the  derivative,  r‘ (0),  If  It  exists.  Is 
necessarily  zero  ( r*  ( 5)  =  -  \  W(k)  Jj  (k5)  k^dk 
=  0,  at  5  =  0):  Thus,  we  are  forced  to  let 
r  (5)  *  1  +  r"(0)52/2  which  gives  an  Incorrect 
prediction  when  compared  with  most  of  the  ex¬ 
perimental  results.  The  maximum  value  of  5 
exp  [  -K (1  -  r  ( 5 ) )]  really  occurs  at  a  point 
50  >0.  It  Is  easy  to  see,  however,  that  50 
becomes  smaller  as  K  gets  larger.  Since  we 
are  now  considering  the  case  when  K  Is  large, 

50  Is  necessarily  a  small  number.  Hence,  the 
proper  approximation  to  r  (5)  Is 

r(5)  =  1  +  r'  (50)  5. 

Note  that  r*  (50)  Is  not  zero  since  r  ( 5)  Is 
monotone  decreasing  at  least  for  small  values 
of  5.  If  so,  the  Integrated  result  of  (1)  Is 

P(0)  -  — -7-  C  (  cos^O  +  C  sln20)  ”3^  2  (4) 

16/ 

where  C  =  (\/  [47ro2r'  (€0)1>2.  It  Is  Important 
to  observe  that  C  does  not  vary  as  x  .  The 
reason  Is  that  5Q  varies  as  K  varies  l.e. , 
as  X,  o.  or  0  varies.  Comparison  with  experi¬ 
mental  results4,  3  Indicates  that  C  should  de¬ 
crease  as  X  decreases.  There  Is,  however,  no 
definite  pattern  for  C  vs  X  since  for  different 
surfaces  r*  (5)  may  be  different.  When  we  ex¬ 
press  r ( 5)  In  terms  of  W(k),  we  see  that  the 
magnitude  of  r*  (5)  Is  larger  If  the  surface  has 
more  high  frequency  components.  If  a  rougher 
surface  Is  defined  to  have  more  high  frequency 
components,  then  a  rougher  surface  calls  for  a 
smaller  C.  We  may  also  want  to  say  that  the 
surface  Is  rougher  If  Its  o  Is  larger.  If  so, 
smaller  C  Indicates  a  rougher  surface  In  both 
cases.  From  the  dependence  of  C  upon  o  and 


r‘  (Cq)»  It  Is  clear  that  this  theory  Is  difficult  to 
pin  down  J>y  experiments  unless  we  know  the 
Surface  in  great  details.  A  general  conclusion 
for  the  case  when  K  Is  large  is  that  P  (0)  tends 
to  be  independent  of  X  and  gets  flatter  with 
angular  change  as  frequency  Increases. 

Comparison  between  the  two  cases  shows 
that  It  Is  easier  to  check  the  surface  roughness 
spectrum  for  the  case  X»o  since  the  effect  of 
a  Is  Independent  of  the  variations  of  W(k'). 
When  K  is  large,  <y2  acts  together  with  r*  (50) 
to  determine  C.  This  may  lead  to  unexpected 
results.  For  example,  we  may  have  oi>a2  and 
OJ  Uol)  <r2*(502)  resulting  In  (£0,)  < 

ff2 *2'£q2) •  Thus  we  get  Ci  >  C2  although  oj 
>  02-  Consequently,  P(0)  curves  for  surface 
no.  1,  which  has  larger  standard  deviation  than 
surface  no.  2,  drop  off  faster  with  Incident  angle. 

The  above  discussion  for  both  cases 
shows  that  a  precise  theoretical  prediction  on 
frequency  behavior  using  Klrchhoff's  theory  is 
difficult.  However,  meaningful  qualitative 
picture  may  be  obtained.  When  the  properties  of 
the  surface  under  study  Is  known,  verification  of 
the  theory  is  possible^  No  verification  was 
possible  in  the  past1,3  since  r‘ (0)  was  used  In 
place  of  r1  (5„)  and  r'  (0)  was  Interpreted  as  one 
over  the  correlation  distance. 

COMPARISON  WITH  EXPERIMENTAL  RESULTS 

Two  model  surfaces  with  different  statis¬ 
tical  properties  are  considered.  Surface  no.  1  Is 
made  of  rubber  polymer  with  a  clay  mold  and 
surface  no.  2  is  made  of  steel.  Both  surfaces 
have  gentle  slopes  so  that  Klrchhoff's  theory  Is 
applicable.  Profiles  of  these  surfaces  are 
measured  with  a  depth  gauge  at  Intervals  0.25  mm 
and  0.5  mm  apart  for  surfaces  no.  2  and  1, 
respectively.  The  standard  deviations,  the  cor¬ 
relation  functions  (see  Fig.  1  and  2)  and  the 
slopes  of  the  correlation  functions  (see  Fig.  5) 
are  then  computed  from  the  profiles  obtained. 

The  acoustic  simulator  facilities  of  the 
Center  for  Research  at  the  University  of  Kansas 
are  employed  for  measuring  the  average  back- 
scattered  power.  This  Is  done  at  frequencies 
ranging  from  0.5  to  1.5  MHz  (see  Fig.  3  and  4). 
Only  Incident  angles  In  the  range  0*  to  30* 
measured  from  the  vertical  are  considered,  since 
there  Is  no  general  agreement  as  to  the  validity 
of  the  theory  for  large  angles  of  Incidence.  (It 
Is  not  clear  how  large  the  Incident  angle  can  be 
before  the  theory  breaks  down).  Also,  for  this 
range  of  Incident  angles  we  are  sure  that  the 
scattering  behavior  of  the  acoustic  waves  is  the 
same  as  that  of  the  electromagnetic  waves. 

Hence,  the  conclusions  arrived  at  In  this  paper 
applies  equally  well  to  electromagnetic  scatter¬ 
ing  . 
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To  check  the  theory  given  by  (4),  we  pro¬ 
ceed  as  follows.  First  of  all.  for  each  frequency 
we  can  fit  (4)  to  the  corresponding  experimental 
curve  to  determine  the  parameter  C.  Secondly, 
we  calculate  r* (Co)  *  x/(4jr cr  C1'2 ).  Thirdly, 
we  determine  ?G  as  the  point  where  the  factor 
{  exp  (-  K  (1  -r(?))]  ,  is  maximum.  The  cosz0 
factor  that  appears  in  K  is  evaluated  at 
0  •*  IS*  assuming  C  gives  the  closest  fit  at 
0  *  15*.  Actually,  the  effect  due  to  cos  0  is 
not  significant  as  compared  with  r  ( f )  at  least 
for  0  between  0*  and  20* .  The  r'  (?0)  thus 
obtained  is  our  predicted  slope  for  r'  (?)  at  ?0. 
Now  on  the  r'  (?)  curve  calculated  from  the  pro¬ 
file  of  the  surface,  we  can  also  find  a  value  for 
r'  (?)  at  ?g.  Let  us  call  this  r' (?D)  our 
measured  slope  as  opposed  to  the  predicted 
slope  for  r(?).  Comparison  of  the  predicted 
with  measured  values  of  r'  (?„)  at  different 
frequencies  is  shown  in  Fig.  5.  Very  good  agree¬ 
ment  Is  obtained.  The  result  for  surface  no.  2  is 
more  satisfactory,  however,  since  the  slopes  are 
larger.  As  for  surface  no.  1,  its  slopes  are  of 
the  size,  10"2  -  10.  Hence,  a  small  error  in 
the  third  place  after  the  decimal  in  evaluating 
the  slopes  from  the  profile-data  will  produce 
tferge/percentage  error.  However,  we  can  still 
see  that  the  order  of  magnitude  is  correct  and 
the  agreement  Is  within  experimental  error. 


CONCLUSION 


Measurements  for  the  backscattered 
power  from  statistically  known  rough  surfaces 
are  made  at  various  frequencies.  Results  are 
used  to  verify  the  scattering  law  obtained  under 
Klrchhoff's  approximation.  When  the  parameters 
In  the  scattering  law  are  properly  Interpreted 
good  agreement  Is  obtained. 


ACKNOWLEDGEMENT 

This  work  is  supported  by  the  National 
Science  Foundation  under  grants  GK  875  and 
GK  1153. 


REFERENCES 


1.  Hagfors,  T.,  "Relationship  of  Geometric 

Optics  and  Autocorrelation  Approaches  to 
the  Analysis  of  Lunar  and  Planetary  Radar," 
T.  Geophys.  Res. .  vol.  71,  p.  379- 
383,  January  1966. 

2.  Rice,  S.  O.,  "Reflection  of  Electromagnetic 

Waves  from  Slightly  Rough  Surfaces,  " 
Comm.  Pure  Appl.  Math  4.  pp.  351- 
378,  1951. 

3.  Beckmann,  P.,  "Scattering  by  Composite 

Rough  Surfaces."  IEEE,  vol.  53, 
pp.  1012-15,  August  1965. 

4.  Chapman,  R.  P.  and  H.  D.  Scott,  "Surface 

Backscattering  Strengths  Measured  Over 
an  Extended  Range  of  Frequencies  and 
Grazing  Angles, "  I.  Acous.  Soc,  Am. . 
vol.  36,  pp.  1735-1737,  September  1964. 

5.  Beckmann,  P.  and  W.  K.  Klemperer,  "Inter¬ 

pretation  of  the  Angular  Dependence  of 
Backscattering  from  the  Moon  and  Venus," 
Radl9.SslgrK9>  vol.  69D, 
pp.  1669-76,  December  1965. 


Figure  1.  Statistical  Properties  of  Surface  No.l  Figure  2.  Statistical  Properties  of  Surface  No. 2 
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Figure  4.  Average  Relative  Backscattered  Power 
vs  Angle  of  Incidence  from  Surface 
No.  2 


Figure  3.  Average  Relative  Backscattered  Power 
vs  Angle  of  Incidence  from  Surface 
No.  I 


Figure  5.  Comparison  Between  Predicted  and 
Measured  Slopes  of  the  Correlation 
Functions  for  Surfaces  No.  1  and  No.  2 
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I.  RESULTS  OF  SCATTERING  EXPERIMENT 

The  scattering  experiment  consisted  of  the  measurement  of  back- 
scatter  from  two  statistically  rough  surfaces.  The  fundamental  surface 
was  gently  undulating  surface  having  slopes  large  relative  to  the  incident 
wavelength.  The  other  surface  was  produced  by  adding  sand  particles  to 
the  fundamental  surface.  As  shown  by  the  results,  the  effect  of  super¬ 
imposing  small  scatterers  onto  the  fundamental  surface  was  to  convert 
the  frequency  dependent  scattering  coefficient  curve  from  a  negative 
slope  to  a  positive  slope.  An  attempt  was  made  to  explain  this  and  other 
behavioral  characteristics  and  to  relate  the  results  to  theoretical  predic¬ 
tions  of  frequency  dependence. 

The  experiment  results  presented  in  this  chapter  were  recorded  as 
described  in  Chapter  4  of  Reference ^  and,  to  the  degree  possible  by 
repetitive  measurements,  are  considered  to  be  accurate  to  within  equipment 
tolerances.  Data  that  could  not  be  repeated  were  not  included  in  the 
results.  These  include  a  sharply  rising  behavior  of  the  scattering  coefficient 
curve  in  the  region  2.5  MHz  to  3.0  MHz  discovered  in  the  measurement  of 
the  fundamental  surface,  which  could  not  be  checked  due  to  subsequent  mal¬ 
function  of  the  transducer  pair  used  in  this  frequency  region.  This  mal¬ 
function,  and  the  unreliable  behavior  of  the  0.3  MHz  to  0.1  MHz  region 
of  the  low-frequency  transducer  pair,  restricted  the  range  of  the  available 
data. 

A.  Measurement  Results  of  Fundamental  Surface 

The  behavior  of  the  scattering  coefficient  as  a  functio  .  of  frequency 
for  the  fundamental,  gently  undulating  rough  surface  is  shown  in  figures 
1  through  5  at  incidence  angles  of  0° ,  15"  ,  30°  ,  45° ,  60°  .  Figure  6  shows 
all  the  curves  plotted  together. 
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The  curves  show  two  behavioral  regions;  one  below  and  one  above 
about  1 .25  MHz  (\  =  1 .2mm) .  In  the  region  below  1  .25  MHz  the  0*  and 
15°  data  have  a  negative  slope  which  increases  above  1.25  MHz.  For 
30®  ,  45®  ,  and  60®  the  lower  frequency  region,  at  least  down  to  0.5  MHz 
(\  =  3.0mm),  the  slope  is  approximately  zero.  Below  0.5  MHz  a  positive 
slope  is  suggested,  especially  at  45®  and  60®,  however  the  confidence 
level  of  these  data  is  not  as  high  as  that  above  0.5  MHz. 

The  region  above  1.25  MHz  has  a  pronounced  negative  slope  for 
all  angles.  This  slope  is  confirmed  to  2.5  MHz  (X  =  0.6mm)  and  is 
believed  tc  continue  to  higher  frequencies,  however  some  measurements 
indicated  a  rapid  change  to  a  large  positive  slope  in  the  region  above 
2.5  MHz.  These  measurements  could  not  be  confirmed  and  no  physical 
explanation  was  found  for  such  behavior,  therefore  the  data  were  dis¬ 
carded  . 

The  data  were  also  plotted  as  a  function  of  angle  for  several  dis¬ 
crete  frequencies  in  figure  7.  The  plots  indicate  a  greater  negative 
slope  near  vertical  for  the  low  frequencies  with  a  gradual  flattening  out 
to  about  1  .2  MHz  (X  =  1  .2mm) .  Thereafter  the  plots  hold  approximately 
the  same  shape  but  gradually  drop  in  overall  amplitude. 

B.  Measurement  Results  of  Sanded  Surface 

The  frequency  dependence  of  the  scattering  coefficient  measured 
for  the  fundamental  surface  densely  coated  with  0 .595-0 .825mm  sand 
particles  are  shown  in  figures  8  through  12  at  incidence  angles  of  0® ,  15®  , 
30° ,  45®  ,  and  60®  .  Figure  13  shows  all  curves  plotted  on  the  same  graph. 

As  in  the  case  of  the  fundamental  surface,  these  data  exhibit  dif¬ 
ferent  behavior  in  two  frequency  regions.  In  the  region  below  1.25  MHz 
the  curves  have  a  positive  slope  of  approximately  the  same  magnitude 
for  all  angles.  Above  1 .25  MHz  the  data  tend  to  become  frequency  inde¬ 
pendent.  There  was  some  support  for  defining  a  third  region  above  about 
2.0  MHz  (X  =  0.75mm)  where  the  data  exhibit  a  positive  slope  such  as 
was  found  for  the  fundamental  surface.  Unfortunately  this  could  not  be 
confirmed  due  to  the  malfunction  of  the  transducer  pair  used  in  this  region. 
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These  data  were  also  plotted  as  a  function  of  angle  for  several 
discrete  frequencies  in  figure  14.  The  plots  indicate  a  greater  nega¬ 
tive  slope  near  vertical  for  the  low  frequencies  with  a  gradual  flattening 
as  the  frequency  increases  out  to  about  1 .25  MHz  (\  =  1 .2mm)  just  as 
in  the  case  of  the  fundamental  surface.  The  measurements  at  45*  inci¬ 
dence  confuse  the  interpretation  of  the  data  and  an  immediate  explanation 
Is  not  available  for  this  phenomenon.  The  readings  tend  to  conform  to 
expectations  as  the  frequency  increases.  Above  1.25  MHz  the  curves 
remain  almost  unchanged  within  the  equipment  accuracy  tolerances. 

C .  Comparison  of  Results  with  Experiments 

(2) 

Following  the  analysis  procedure  used  extensively  by  Katz'  ' 
the  scattering  coefficient  frequency  dependence  curves  were  fit  to  a  func¬ 
tion:  a*  =  a\°.  The  fit  was  done  by  frequency  regions  and  a  determination 
of  the  best  o  was  made  for  all  data.  In  figures  15  and  16  the  variation  of 
the  coefficient  o  is  plotted  as  a  function  of  incidence  angle  for  the  funda¬ 
mental  surface.  In  figures  17  and  18  the  coefficients  are  plotted  for  the 
sanded  surfaces. 

In  the  next  section  an  attempt  will  be  made  to  explain  the  behavior 

of  the  curves  using  the  guidelines  offered  by  the  Kirchhoff  theoretical 

method  of  expressing  the  scattering  processes,  however  as  reference  to 

other  experimental  results,  the  values  of  the  coefficient  a  for  angular 

backscatter  measurements  of  sea  clutter,  cities,  snow  covered  terrain, 

grass  surfaces,  the  moon,  and  asphalt  and  concrete  roads  are  shown  In 

(2) 

figure  19.  These  data  were  compiled  by  Katz'  '  using  data  from  the  Ohio 
State  University  and  Naval  Research  Laboratory  measurements  programs. 
These  data  were  obtained  by  polychromatic  systems,  hence  the  slope  factors 
were  obtained  using  widely  spaced  points.  Such  subtleties  as  were  found  in 
the  measurements  of  the  surfaces  used  in  this  program  could  not  have  been 
detected . 

D .  Comparison  with  Theory 

Although  some  degree  of  matching  of  the  measured  and  theoretical 
results  was  obtained,  considerable  doubt  exists  as  to  the  ability  of  the 
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Kirchhoff  method  to  predict  the  frequency  dependence  of  backscatter,  at 

least  in  cases  where  the  actual  surface  autocorrelation  is  approximated  by 

an  autocorrelation  function.  A  point  of  major  concern  is  that  the  best  fit 

occurs  using  the  expressions  derived  for  an  exponential  autocorrelation 

function,  which  clearly  cannot  fit  the  measured  surface  autocorrelation. 

This  contradiction  is  somewhat  overcome  by  the  approach  used  by  Fung  and 
(3) 

Leovaris  ' ,  which  allows  incorporation  of  autocorrelation  functions  other 

than  the  exponential.  Even  with  this  approach  the  maximum  wavelength 
+2 

variation  is  X  (at  ©  =  0® ) ,  and  the  coefficient  must  decrease  as  ©  in¬ 
creases.  This  characteristic  is  contradicted  by  the  measurements  of  the 

(4) 

fundamental  surface  (see  figure  6)  and  the  results  of  Wiltse,  et  al.  '. 

Likewise  the  theory  cannot  predict  results  such  as  those  compiled  by  Katz. 

The  effective  parameter  approach  is  supported  by  the  physical 

reasoning  that  different  subranges  of  structure  sizes  dominate  the  return  at 

different  incidence  angles  and/or  wavelengths  .  These  subranges  have 

autocorrelations  and  variances  which  are  not  the  same  as  the  autocorrelation 

and  variance  of  all  subranges  combined,  i.e.  the  entire  surface.  Several 

attempts  have  been  made  to  incorporate  this  character  into  the  scattering 

(5  6  7  8) 

expressions  to  describe  the  angle  dependence'  '  '  '  .  These  attempts 

have  been  only  moderately  successful. 

In  summary,  the  measurements  obtained  in  this  program,  as  well  as 

(8) 

those  of  other  programs,  suggest,  in  agreement  with  Davies'  ,  that 
use  of  the  autocorrelation  function  approximation  of  the  surface  autocorrela¬ 
tion  does  not  give  sufficient  statistical  information  to  describe  the  scattering 
from  the  surface. 

II.  IMAGE  EXPERIMENT 

In  recent  years  the  use  of  imaging  radar  systems  for  both  military 
reconnaissance  and  geoscience  investigations  has  expanded  rapidly. 
Associated  with  this  expanding  demand  has  been  the  growing  realization 
that  the  monochromatic  nature  of  such  systems  handicaps  their  utility, 
especially  in  military  reconnaissance  applications.  This  chapter  presents 
an  experiment  with  a  panchromatic  imaging  system.  In  this  experiment 
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monochromatic  images  of  complex  targets  are  compared  with  panchromatic 
images  to  show  the  effect  of  frequency  averaging.  The  substantial  improve¬ 
ment  in  target  definition  obtained  is  impressive,  especially  in  view  of  the 
fact  that  less  than  a  +  10%  bandwidth  is  employed,  and  that  little  increase 
in  system  power  or  complexity  is  required  to  gain  this  improvement. 

The  theme  of  the  chapter  is  illustrated  in  Figure  20  which  shows 
an  acoustic  image  of  two  spheres.  The  images  actually  show  an  expanded 
view  of  a  small  area  of  the  sphere's  surface  near  the  apex.  The  return 
signal  strength  from  regions  away  from  the  apex  is  below  the  noise  level. 

In  20  (a)  the  spheres  are  imaged  monochromatically  and  lobing  is  evident. 

In  20  (b)  the  frequency  is  swept  +10%  from  the  center  frequency.  The 
frequency  averaging  of  the  lobes  is  obvious.  The  monochromatic  lobes 
are  investigated  in  detail  in  section  II-C  and  the  effect  of  frequency  aver¬ 
aging  is  illustrited. 

A.  Measurement  Technique 

Acoustic  images  were  obtained  in  a  manner  completely  analogous 
to  an  airborne  imaging  radar  system.  The  antennas  (transmitting  and 
receiving  transducers)  were  moved  past  the  targets  at  a  fixed  distance 
and  uniform  velocity.  The  transducers  were  mounted  on  a  linear  motion 
carriage  that  was  capable  of  a  wide  range  of  movement  velocities  and 
transducer  position  settings.  The  position  geometry  is  shown  in  Figure  21  and 
a  photograph  of  the  arrangement  is  also  shown. 

The  system  operation  associated  with  obtaining  acoustic  images 
was  as  follows:  The  signal  return  is  envelope  detected  and  applied  to  the 
grid  of  an  imaging  CRT.  Each  return,  i.e.  the  return  associated  with  each 
transmitted  pulse,  intensity  modulates  one  horizontal  line  trace  on  the  CRT. 
Sequential  returns  form  sequential  intensity-modulated  line  traces  each  vertically 
displaced  from  each  other.  The  end  result  is  an  intensity  modulated  raster 
display.  A  photographic  image  is  obtained  by  exposing  the  film  during  the 
entire  vertical  traverse  time  interval.  The  producure  is  illustrated  in 
figure  22 . 

The  imaging  experiment  was  designed  around  the  limitation  presented 
by  the  transducers.  That  is,  the  available  piston  transducers  provide  only 
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conical  beams  and  have  relatively  high  sidelobes.  The  use  of  spherical 
targets  insured  that  the  sidelobes  of  the  antenna  did  not  influence  the  return 
over  the  area  of  investigation  of  the  reradiation  pattern.  Ideally  the  experi¬ 
ment  would  concentrate  on  imaging  area -extensive  targets  such  as  those 
of  concern  to  imaging  radar  systems,  however  the  lack  of  a  fan-beam  antenna 
discouraged  this  approach.  Fortunately,  as  will  be  shown  in  section  II-B, 
the  spherical  targets  provided  the  information  needed  to  relate  to  the  area 
extensive  imaging  problem,  hence  the  limitation  to  a  conical  beam  did  not 
restrict  the  value  of  the  experiment. 

The  parameters  of  the  experiment  were  set  with  the  following  con¬ 
sideration.  The  objective  of  the  experiment  was  to  determine  the  effect 
of  introducing  a  broad -spectrum  signal  in  place  of  a  monochromatic  signal. 

The  targets  examined  were  complex  and  consequently  provide  a  reradiation 
pattern  consisting  of  a  multiple  of  fine  lobes.  To  determine  the  effect  of 
frequency  averaging  required  that  the  effect  of  sample  averaging  be  excluded. 
That  is,  the  image  would  be  altered  if  the  number  of  "looks"  were  changed 
by  the  fact  that  more  or  less  lobes  would  be  seen.  To  elimir^e  this  variable, 
the  sampling  rate  (PRF)  was  set  sufficiently  high  that  all  lobes  would  be 
recorded  and  the  rate  was  not  altered  when  the  swept  frequency  was  intro¬ 
duced.  The  darriage  velocity  was  set  at  6.5  mm/sec  and  the  PRF  was  set 
at  approximately  130  Hz.  This  provided  one  sample  for  each  .05mm  of 
carriage  travel. 

The  system  employs  the  "slow-sweep"  mode.  This  mode  of  operation 
consists  of  sampling  the  swept  frequency  at  the  pulse-rate  frequency.  Ideally 
broad -spectrum  imaging  would  be  done  using  a  "fast-sweep"  mode,  that  is 
the  entire  spectrum  would  be  transmitted  during  each  pulse.  This  could 
conceivably  be  simulated  by  increasing  the  PRF  proportional  to  the  sweep 
time.  That  is,  if  the  PRF  for  monochromatic  imaging  os  100  Hz,  set  the 
sweep  time  at  (1/100)  sec  and  increase  the  PRF  to  400  Hz  to  provide  4  samples 
per  sweep.  Unfortunately  this  reintroduces  the  effect  of  increasing  the 
sampling  rate  and  hence  for  this  experiment  the  PRF  for  both  monochromatic 
and  swept  frequency  signals  was  the  same,  130  Hz.  The  sweep  time  was 
adjusted  to  give  4  samples  per  sweep  (about  23.5  m  sec).  The  samples 


B7 


were  not  synchronized  with  the  sweep,  and  hence  the  frequency  sampling 
varies  sweep-to-sweep.  In  this  manner  the  full  sweep  is  effective  even 
though  the  sampling  rate  is  very  low. 

The  width  of  the  pulse  was  set  at  .6  -  .8  m  sec  to  provide  a  return 
segment  which  did  not  include  the  effect  of  the  high  frequency  components 
In  the  leading  and  trailing  edge  of  the  transmitted  pulse.  It  was  determined 
that  the  spherical  targets  caused  the  return  signal  to  be  restricted  to  an 
area  near  the  apex,  hence  the  leading  edge  effect  decayed  rapidly.  There¬ 
fore  by  extending  the  pulse  width  beyond  this  decay  interval,  the  leading 
edge  effect  was  excluded .  This  effect  is  evident  in  the  images  shown  in 
figure  20 . 

The  useful  bandwidth  of  the  swept  signal  was  restricted  by  the 
frequency  response  of  the  transducers.  The  response  curves  fall  so  rapidly 
that  variations  greater  than  about  +  10%  from  resonance  were  impractical. 

B.  Targets 

Three  primary  targets  were  used  in  this  experiment:  (1)  a  4.5  cm 
diameter  rubber  sphere;  (2)  a  7.8  cm  diameter  styrofoam  sphere;  and  (3) 
an  8-element  planar  array  consisting  of  2 .5  cm  diameter  wooden  spheres 
located  on  3.2  cm  centers  in  two  rows  of  4-elements  each.  The  targets 
are  shown  in  the  photograph  as  they  were  positioned  in  the  tank.  The 
rubber  sphere  is  on  the  left. 

1 .  Rubber  Sphere  .  This  target  is  a  relatively  smooth,  soft  rubber 

sphere  selected  to  provide  a  broad,  hopefully  frequency-insensitive,  lobe 

in  the  direction  of  backscatter.  Although  the  exact  composition  of  the  sphere 

is  unknown,  it  is  reasonable  to  assume  negligible  shear  waves  and  a  reasonably 

high  reflection  coefficient.  The  diameter  of  the  sphere  is  very  large  relative 

to  the  incident  wavelength,  but  is  smaller  than  the  illuminated  area  of  the 

beam.  Under  these  conditions  the  scattering  cross  section  should  be  almost 

(9) 

independent  of  wavelength  and  the  sphere  would  provide  the  necessary 
calibration  reference  for  other  measurements. 

The  measurements  presented  in  section  II-C  verify  the  presence  of 
a  strong,  almost  frequency  insensitive,  reflection  in  the  backscatter  direc- 
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ticn,  however  the  fact  that  the  sphere  is  not  perfectly  smooth  to  the  incident 
wavelength  is  evident. 

2.  Styrofoam  Sphere.  This  target  was  selected  to  illustrate  an  effect 
very  common  to  imaging  systems  but  one  that  is  difficult  to  simulate  acoustically 
in  a  controlled  manner;  that  is,  an  area  extensive  surface  consisting  of 
various  size  randomly  oriented  scatterers.  The  surface  of  the  sphere  is  rough 

in  the  somewhat  facet-like  manner  typical  of  styrofoam.  Consequently  some 
scatterers  are  sufficiently  close  together  to  form  highly  directional  reradi- 
ation  patterns  similar  to  closely  spaced  arrays.  In  addition  the  multilobe 
pattern  of  widely  spaced  elements  exists,  and  scatterer  sizes  both  large 
and  small  relative  to  the  incident  wavelength  are  present.  The  spherical 
configuration  reduces  the  illumination  area  edge  effects  and  helps  accent 
the  behavior  of  the  surface  reradiation  lobes  by  restricting  the  contributing 
area  to  a  much  smaller  size  than  is  otherwise  possible  using  the  available 
acoustic  transducers  . 

3 .  Array.  The  array  target  consists  of  eight  spheres  arranged  as 
shown  in  figure  23.  The  purpose  of  this  target  is  to  provide  a  large  number 
of  narrow  lobes  of  different  amplitudes.  The  array  is  an  extension  of  a 
4-element  array  (one  row  of  spheres)  which  provided  the  desired  fine 

lobe  structure,  but  the  lobes  were  of  almost  equal  amplitude.  The  effect 
of  variations  in  incident  wavelength  is  better  illustrated  by  use  of  the 
more  directional  pattern  of  the  8-element  array. 

The  reradiation  pattern  of  an  array  of  point  scatterers  spaced  several 
wavelengths  apart  can  be  calculated,  and  the  results  show  a  number  of  lobes 
proportional  to  the  number  of  elements .  The  width  of  each  lobe  is  also  pro¬ 
portional  to  the  number  of  elements .  The  spherical  elements  of  the  target 
array  cannot  be  interpreted  as  point  scatterers;  however  the  area  on  the  sphere 
contributing  to  the  reradiation  is  sufficiently  small  that  a  very  similar  pattern 
to  that  calculated  is  created. 

C .  Results 

The  results  clearly  indicate  that  the  addition  of  frequency  averaging 
of  the  reradiation  pattern  of  complex  targets  substantially  improves  the  de- 
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finition  of  the  target.  The  reader  should  bear  in  mind  that  the  results  pre¬ 
sented  represent  only  the  addition  of  a  swept  frequency  during  transmission, 
done  in  such  a  manner  that  each  pulse  modulated  a  different  carrier  fre¬ 
quency  on  a  repetitive  sweep  basis.  Several  other  techniques  are  possible, 

(11  12  131 

some  perhaps  more  optimum  than  the  method  used  here  '  '  ' .  However, 

the  point  is  well  made  that  the  addition  of  frequency  averaging  warrants  con¬ 
sideration  in  the  design  of  imaging  radar  systems. 

1.  Rubber  Sphere .  Figure  24  shows  the  acoustic  image  of  the  rubber 
sphere  for  four  separate  frequencies  and  a  swept  frequency.  The  images 
were  recorded  from  a  distance  of  51  cm.,  and  consequently  the  illumina¬ 
tion  region  of  the  transmitted  beam  was  6  cm  in  diameter.  The  target  dia¬ 
meter  was  4.5  cm. 

The  images  of  the  styrofoam  sphere  and  the  array  show  the  exist¬ 
ence  of,  what  may  properly  be  termed,  multiple  lobes.  The  return  from 
the  rubber  sphere  consists  of  multiple  "specular  points."  That  is,  the 
sphere  is  obviously  not  perfectly  smooth  to  the  incident  wavelength, 
but  the  roughness  does  not  produce  an  array  effect  such  as  the  styrofoam 
sphere.  The  significance  of  this  distinction  is  that  for  an  array  the  lobes 
“move"  with  frequency  and  hence  an  averaging  is  obtained.  This  effect 
does  not  occur  with  the  rubber  sphere.  The  distinction  evident  between 
images  at  different  single  frequencies  is  that  each  "specular  point"  has 
a  different  amplitude  from  frequency-to-frequency .  The  swept  frequency 
image  clearly  shows  the  target  to  have  two  major  "specular  points",  where 
one  is  6  db  greater  than  the  other.  Other  such  points  are  indicated  to 
exist  at  lower  levels.  No  single  frequency  image  clearly  defines  both 
points,  and  limited  to  any  one  single  frequency  image  the  definition  of 
the  target  is  misleading.  In  comparing  figure  24  to  figure  20,  both  of  which 
show  the  rubber  sphere,  the  reader  should  note  that  figure  20  was  obtained 
by  imaging  a  region  of  the  sphere  slightly  off-center  and  hence  the  dynamic 
range  of  the  "specular  point"  return  was  less  than  the  returns  shown  in  figure 
24,  which  are  from  the  apex. 


2.  Styrofoam  Sphere.  Figure  2  5  shows  the  acoustic  image  of  the 
styrofoam  sphere  for  several  single  frequencies.  In  figure  26  one  of  the  mono- 
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chromatic  images  is  compared  to  swept  frequency  images  of  +_  3%  and  +_  6% 
from  center  frequency.  The  images  were  recorded  at  a  distance  of  approxi¬ 
mately  1  m  and  a  beam  illumination  region  of  11.6  cm  diameter.  The  re¬ 
radiation  lobe  pattern  is  more  complex  than  the  simple  rubber  sphere  and 
the  effect  of  varying  frequency  is  marked.  Of  special  interest  is  the  phase 
cancellation  of  the  center  main  lobe  (well  defined  at  1520  KHz  and  1626  KHz) 
at  1427  KHz,  1744  KHz,  and  1816  KHz.  The  effect  occurred  at  numerous 
intermediate  frequencies  as  well. 

The  application  of  swept  frequency  illumination  (figure  26)  smoothed 
the  monochromatic  lobes  and  formed  a  uniform,  broad  backscatter  lobe. 

The  effect  of  increasing  the  bandwidth  was  to  strengthen  the  off-center  lobes 
found  at  the  higher  and  lower  frequencies  and  hence  "fill"  the  sphere  return. 
Unfortunately  at  sweeps  abov*-  about +_  10%  the  frequency  response  curve 
of  the  transducers  has  decreased  such  that  only  a  small  contribution  is 
obtained  from  the  extreme  frequencies  of  the  sweep,  hence  the  full  effect 
of  frequency  averaging  could  not  be  investigated.  This  is  especially  evident 
in  the  array  target  data . 

3.  Array.  Figure  27  shows  the  acoustic  images  of  the  8-element  array 
at  three  monochromatic  frequencies.  The  images  were  recorded  for  several 
gain  settings  to  show  the  various  amplitude  levels  of  the  multiple  lobes. 

As  the  gain  is  increased  the  lower  level  lobes  become  evident.  Because 
of  the  gain  level  sensitivity  and  the  large  number  of  lobes  it  is  difficult 
to  determine  the  difference  in  the  lobe  structure  at  different  frequencies. 

Figure  28  compares  three  monochromatic  images  at  approximately  the  same 
loop-gain  levels,  i.e.  the  receiver  gain  is  set  to  compensate  for  the  de¬ 
crease  in  transducer  gain. 

Figure  29  shows  the  effect  of  adding  a  swept  frequency  signal 
of  two  different  sweep  widths:  +  100  KHz  and  +_  50  KHz.  As  in  figure  27, 
the  lobe  levels  were  interrogated  by  varying  the  receiver  gain.  The  +  50 
KHz  swept  signal  was  inadequate  to  completely  smooth  the  reradiation 
pattern,  however  the  +100  KHz  swept  signal  provided  good  smoothing  at 
nearly  all  intensity  levels.  It  is  evident  that  a  greater  sweep  width  is 
required  to  completely  "fill"  the  image. 
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Figure  7  Variation  of  scattering  coefficient  with  angle  and 
frequency  -  fundamental  surface. 
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Figure  17.  Variation  of  slope  for  sanded  surface  data 
below  1.25  MHz 


Figure  18.  Variation  of  slope  for  sanded  surface  data 
above  1.25  MHz 
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(a)  Wavelength  dependence  of  grass, 
concrete,  and  asphalt  surfaces. 
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(b)  Wavelength  dependence  of  sea 
clutter.  This  figure  covers  the 
microwave  region  from  0.86  to 
71  cm  and  applies  to  rough  sea 
conditions . 
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(c)  Wavelength  dependence  of  the 
moon's  surface. 


Figure  19 


Variation  of  slopes  for  natural  surfaces 
(after  Katz  1966) 
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Effe  t  ;  receiver  gam  variations  on  images  of  8-element 
arras’  at  three  monochromatic  frequencies 
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Effect  f  swept  frequency  on  image  of  6-olement  array 
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SYSTEM  AND  COMPONENT  CONSIDERATIONS  FOR  THE 
POLYPANCHROMATIC  RADAR 


P.  D.  Shew 


Technical  Memorandum  133-1-2 

SYSTEM  AND  COMPONENT  CONSIDERATIONS  FOR  THE 
POLYPANCHROMATIC  RADAR 

by 

P.  D.  Shaw 


1.0  INTRODUCTION 

The  polypanchroma  tic  radar  is  a  dovlco  for  performing  bioad  spectrum 
electromagnetic  scattering  and  Imaging  measurements  .  The  calculations 
necessary  to  define  the  transmitter  power  and  antenna  wore  previously 
performed  (Reference  1).  Since  publication  of  this  memorandum  the  systom 
functional  arrangement  has  beon  changed.  Consideration  will  be  given 
to  the  revised  system  and  the  components  necessary  to  lmplomert  It.  Tho 
post  detection  processing  (Reference  3)  and  the  digital  control  system 
(Reference  4)  are  described  separately. 

2.0  SYSTEM  REVIEW 

The  basic  design  goals  have  been  described  (Reference  1).  These 

are: 

1 .  Operation  over  a  bandwidth  of  one  octave  (4-8  GHz) 

2.  Measurement  of  backs catterlng  cross  section  as  a  function  of 

a .  Frequency 

b.  Bandwidth 

c.  Incidence  angle 

3.  Production  of  panchromatic  Imagery  as  a  function  of 

a .  Frequency 

b.  Bandwidth 

These  objectives  are  most  economically  achieved  with  a  slow  sweep, 
pulse  modulated  system  In  which  the  frequency  does  not  change  appreciably 
within  the  pulse  period .  The  frequency  range  Is  divided  into  separated 
narrow  bands.  Full  spectrum  coverage  is  achieved  by  the  overlapping  spectra 
of  the  transmitted  pulses  (Figure  l2)„ 
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Figure  1  shows  tho  functional  units  in  tho  radar  whilo  Figure  2 
shows  tho  componont  arrangoment.  It  will  be  notod  that  tho  components 
havo  boon  dividod  into  two  soctions .  Two  mothoJs  of  construction  aro 
possible . 

1  .  All  components  mountod  on  the  ground  with  the  signal  trans¬ 

mitted  to  the  antenna,  at  the  end  of  tho  boom,  by  coaxial 
cablo  (Rcforonco  1).  This  has  the  advantage  of  componont 
accessibility  and  simplified  duplex  switching,  but  glvos  a 
high  signal  loss  due  to  tho  cablo. 

2  .  The  signal  loss  may  bo  reduced  by  mounting  the  transmitter  and 

receiver  T.W.T.'s  and  associated  components  next  to  tho 
antenna . 

Tho  first  method  of  mounting  will  be  used  initially  to  simplify 
adjustment.  The  second  method  can  be  used  later  to  increase  the  range 
of  the  radar. 

3.0  THE  TRANSMITTER 
3.1  The  T.W.T.  Amplifier 

As  described  (Reference  1)  a  transmitted  power  of  20  watts  is 
required.  The  frequency  of  the  transmitted  signal  must  vary  in  frequency 
from  4  GHz  to  8  GHz.  The  only  suitable  amplifier  is  a  traveling  wave 
tube.  A  T.W.T.  manufactured  by  Alto  Scientific  was  purchased.  Its 
relevant  characteristics  are  shown  (Spec.  1). 

When  using  this  amplifier  it  is  necessary  to  consider  the  harmonic 
power  generated  by  the  tube.  The  most  Important  component  will  be  the 
second  harmonic  generated  when  the  fundamental  is  a  low  frequency.  At 
this  point  both  fundamental  and  second  harmonic  are  within  the  pass  band 
of  the  tube  and  the  second  harmonic  may  be  only  3  to  6  db  less  than  the 
fundamental.  It  is  most  Important  that  these  higher  harmonics  not  be 
reflected  back  into  the  tube  output.  We  must,  therefore,  use  an  Isolator 
at  the  tube  output.  If  these  higher  harmonics  are  to  be  Isolated  from  the 
load  it  will  be  necessary  to  use  a  switched  filter  (Figure  3) . 
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Wo  can  calculato  tho  signal  power  required  from  tho  source  as 
fallows: 

TWT  Gain  ■  58  db  (min.) 

TWT  output  ■  44  dbm  (max.) 

Isolator  insertion  loss  ■  0.4  db  (max.) 

2x  Diodo  Switch  insertion  loss  -  3  db  (max.) 

Signal  Powor  ■  *10.6  dbm 

Ihls  amount  of  power  is  obtained  from  the  source  output  through  a 
directional  coupler. 

3.2  The  Sweep  Frequency  Source 

A  source  is  required  which  is  able  to  rapidly  swcup  the  4*8  GHz  frequen¬ 
cy  range.  The  source  chosen  was  a  swoop  frequency  oscillator  manufac¬ 
tured  by  Hewlett  Packard.  The  relevant  characteristics  are  shown  (Spec.  2) . 

This  oscillator  has  an  internal  P.l.N.  leveling  system  (Section  6.0).  Tho 
relevant  components  were  removed  and  used  externally  to  facilitate  power 
division. 

4.0  THE  RECEIVER 

4.1  Ihg  Pfegmpllftef 

The  requirements  for  the  amplifier  are 

a.  Low  noise  figure  (<  10  db) 

b.  Frequency  range  4-8  GHz 

c.  High  gain 

d.  Reasonable  cost 

A  T.W.T.  Amplifier  is  necessary.  The  specifications  of  the  model 
chosen  are  given  in  (Spec.  3) . 

4.2  Intermediate  Frequency  Amplifier 

An  amplifier  manufactured  by  R.H.G.  Labs  was  chosen  (Spec.  4). 

The  factors  influencing  this  choice  and  the  associated  components  are  des¬ 
cribed. 


C4 


4.2.1  Choice  of  bandwidth 

The  maximum  spectral  width  of  transmitted  signal  Is  obtained  when 
the  pulse  has  the  minimum  width  of  20nS.  This  would  have  a  bandwidth. 

B  -  ^  -  60  MHz 

T 

This  gives  a  sufficient  spectral  overlap  between  the  frequency  shlftod  pulses 
to  approximate  continuous  frequency  coverage. 

The  maximum  range  we  will  assume  to  be  915  m  which  gives  a  signal 
delay  of  6.1  pS .  We  will  allow  100  ns  for  delay  through  the  systom  giving 
a  total  delay  of  6.2ps.  if  the  frequency  source  Is  sweeping  at  Its  maximum 
rate  of  4  GHz  In  10  ms,  the  frequency  will  change  by  2.4  MHz.  We 
will  assume  that  the  doppler  effect  is  negligible  In  our  static  system. 

Tho  bandwidth  of  the  receiver  must  bo  great  enough  to  accommodate  both  the 
bandwidth  of  the  Individual  pulse  and  the  .'requcncy  change  within  tho  delay 
tlmo.  A  receiver  bandwdith  of  100  MHz  was  chosen,  which  will  accomplish 
this  and  still  leave  room  for  possible  future  reduction  of  the  Individual  pulse 
width. 

4.2.2  I.F,  Amplifier  Gain 

We  will  estimate  the  maximum  gain  required  in  the  IF  stage.  Dynamic 
range  will  be  considered  later. 

The  minimum  received  signal 

Pr  «-  KTBFS 

B  =  100  MHz  ■  receiver  bandwidth 
F  =  10  db  =  receiver  noise  figure 
S  =  3  db  =  signal  to  noise  ratio 

If  the  required  I?  power  out  is  1  mw 


Receiver  Gain  =  90.1  db 
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Actual  Receiver  gain  »  32.2  +  23  ♦  55-110.2  db. 
We  must  attenuate  the  received  signal  by  20  db. 


4.3  Receiver  Dynamic  Ranoo 
4.3.1  R.F.  Section 


We  must  first  cstlmato  the  maximum  signal  power  that  appears  at  tho 
receiver  TWT  Input.  This  occurs  when  we  look  at  a  target  at  minimum  range 
end  vertical  incidence. 

We  have 


P  p  2  2  A2 

P.  Pt*  *  d  qo 
r  R2  256 


(Ref.  1) 


Where  Pf  -  return  signal  power  (max.) 

■  20  watts  -  transmitted  power  (max.) 
e  -  0.5  -  antonna  efficiency 
d  -  0.61  m  -  diameter  of  antenna 

aQ  -  20  db  -  radar  scattering  cross  sectlon/unlt  area  (max.) 
R  -11.1m  -range  (rtln.) 


Pf  -  58  mw  -  17.6  dbm 

The  maximum  power  out  of  the  I.F.  preamp  Is  -10  dbm.  to  avoid 
saturating  the  preamp  output.  The  corresponding  gains  are  preamp  23  db. 
and  TWT  32.2  db. 


Max  signal  Into  receiver  TWT  -  -23  -10  -32.2  -  -65.2  dbm. 

We  must  attenuate  the  received  signal  by  17.*  ♦  65.2  -  81 .8  db. 

The  minimum  received  signal  Is  -93.2  dbm 


Dynamic  power  range  of  the  receiver  -  93.2  -  65.2  -  28  db 
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Tho  dynamic  range  sots  the  limits  oi  revived  signal  variation  at 
a  constant  attenuator  sotting.  The  effective  dynamic  range  will  bo  roduced 
since  wo  will  wish  to  uso  finite  steps  of  attonuation.  Theso  rocoived  signal 
‘'ariatlono  are  duo  to  Raylolgh  signal  fading  (Ref.  1)  and  system  parameter 
variations.  We  will  consldor  the  effort  of  those  parameters. 

At  angles  of  incidence  away  from  the  vertical  we  have 


d  2  2  .2 

V  ’  d  ''0 


64  R  sin  €> 


l  Boam  width  limited  mode  ] 


Pt  »  e  d  f  t  o0 

P  - - s -  (  Pulse  length  limited  r.ouo  J 

128  R  sine 


Where  Pf 

pt 

€ 

d 


T 

*o 

f 

e 


return  signal  power 
transmitted  powor 
antenna  efficiency 
diameter  of  antenna 
pulse  length 

radar  scattering  cross  soctlon/unit  area 
signal  frequency 
angle  of  Incidence 


To  obtain  the  scatterometrv  data  we  will  operate  in  tho  beamwidth  limited 
mode  with  approximately  a  constant  range.  The  only  lrequency  dependent 
term  in  the  expression  for  the  return  power  is  the  radar  scattering  cross** 
section.  Since  the  variation  of  this  scattering  cross -section  with  frequency 
is  expected  to  be  small  we  can  obtain  a  set  of  measurements,  at  a  constant 
angle  of  incidence,  whilst  frequency  is  varied,  without  changing  the  R.F. 
attenuator  setting . 

For  the  imaging  measurements  we  will  operate  in  a  pulse  length 
limited  mode  to  improve  the  range  resolution  at  far  range.  For  targets 
where  the  fading  effect  is  great,  the  dynamic  range  may  be  insufficient. 
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It  will  then  be  necessary  to  make  measurements  over  restricted  frequency 
and  range  bands  and  to  superimpose  the  results. 

An  R.F.  attenuator  which  has  steps  of  5  db  and  a  maximum  of  90  db 
will  probably  be  sufficient.  This  can  be  achieved  using  attenuators  of 
5,  10,  20,  40,  30  db  in  series  combinations. 

A  linear  IF  amplifier  does  not  allow  great  dynamic  range .  The  dynamic 
range  can  be  increased,  by  using  a  logarithmic  IF  amplifier.  We  would 
then  be  limited  by  the  noise  power  of  the  receiver  TWT  and  the  maximum 
power  that  can  be  applied  to  its  mixer  load.  It  would  increase  the  dynamic 
range  in  this  case  by  23  db. 

4.3.2  JF  Section 

The  maximum  output  of  the  preamp  is  -lOdbrn  while  that  of  the  post 
amp  is  0  dbm  and  the  post  amp  gain  fs  55  db.  Under  the  condition  of  maxi¬ 
mum  signal  return  we  must  insert  45  db  attenuation  between  the  IF  preamp 
and  IF  postamp.  A  switched  IF  attenuator  has  been  purchased  for  this  pur¬ 
pose. 

4.3.3  Automatic  Gain  Control 

The  automatic  gain  control  input  of  the  amplifier  can  be  used  to 
change  the  amplifier  gain  with  signal  frequency  to  match  the  change  in 
antenna  gain  with  frequency.  This  reduces  the  dynamic  range  require¬ 
ments  of  the  post  detection  components . 

5.0  DUPLEXING 

5.1  Isolation  Requirements 

The  arrangement  of  the  required  switches  and  circulator  are  shown 
in  Figure  4.  The  specifications  of  the  required  components  are  given 
(spec.  5),  (Figures  10,  11). 

Aj  A2  A^  is  the  isolation  of  each  switch  in  db  when  the  switch  is  OFF 
(power  reflected).  A^  is  the  isolation  between  reverse  ports  of  the  circulator. 
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B1  B2  B4  is  insertion  l°ss  when  the  switch  is  ON  (power  transmitted). 
Bg  is  the  insertion  loss  between  forward  ports  of  the  circulator. 

There  are  two  system  states: 

Transmit  S2  ON  OFF 

Receive  S1  Sg  OFF  ON 

On  transmit  it  is  necessary  to  prevent  the  receiver  from  saturating. 
On  receive  it  is  necessary  that  the  power  at  the  receiver  input  due 
to  transmitter  TWT  noise  and  signal  leakage  be  much  less  than  the  received 
radar  return  power. 

We  will  now  define  the  system  parameters: 

Transmitter  TWT  maximum  pwer  out 
Transmitter  TWT  gain 
Transmitter  TWT  noise  figure 
Receiver  TWT  gain 
Receiver  TWT  noise  figure 
IF  Preamplifier  gain 
IF  Preamplifier  power  out 
IF  Post  amplifier  gain 
IF  Post  amplifier  power  out 

On  transmit: 

PIF  max  =  PT  max  +  GR  +  GIF  ”A3  'A4 
a3  +  a4  =  100.2  db 
But  A3  =  18  db 
.*.  A4  =  82.2  db 

But  P.I.N.  diode  switch  isolation  =  2  x  38  db  =  76  db 
We  will  use  10  db  of  fixed  attenuation  between  the 
receiver  TWT  amplifier  and  the  IF  preamplifier. 


P^  max  =  139db* 

Gx  =  43.5  db 

NpT  =  35  db 

GR  =  32.2  db 

NpR  =  7.3  db 

Gjp  =  23  db 

P™  m  =  94  db 
IF  max 

G'Ip  =  55  db 

P'  =  94  db 

IF  max 


*  All  power  is  referenced  to  the  Thermal  noise,  KTB,  where  B  is  the 
receiver  bandwidth  =  100  MHz, 
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To  prevent  the  post  amplifier  from  saturating  we  must 
put  a  diode  switch  between  the  pre  amplifier  and  post 
amplifier  with  an  isolation  of  greater  than  55  db  and  a 
switching  speed  of  less  than  20  nS. 

On  receive: 

We  will  make  the  TWT  noise  and  signal  leakage  power  at 
the  receiver  TWT  input  3  db  above  the  thermal  noise  level. 
Power  due  to  transmitter  TWT  noise  =  NTp  +  GT  -  A2  -  A„ 

A2  +  A3  =  35  +  43.5  -  3  =  75.5  db 

Power  due  to  signal  leakage  =  Prp  mQX  “  “  A2  ~  A3 

Aj  +  A2  +  A3  =  139  -  3  =  136  db 

5.2  Switch  Power  Rating 

The  pulse  to  be  trimsmitted  has  a  pea^  power  to  20  w  and  is  applied 
to  the  switch  S2  only  when  the  switch  is  in  the  open  condition.  The  switch 
S2  must  be  able  to  pass  this  signal  safely  and  without  signal  degradation. 
The  limiting  factor  will  be  diode  voltage  breakdown,  while  the  switch  is  in 
the  open  position,  and  average  power  dissipation.  From  the  specifications 
of  the  switch  we  see  that  the  maximum  insertion  loss  with  zero  bias  is 
1.5  db.  For  this  value  the  average  power  must  not  exceed  10  watts  and  a 
peak  power  of  80  watts. 

Due  to  the  high  cost  of  the  P.I.N.  diode  switches  and  the  many 
possible  failure  modes  we  must  find  some  method  of  protecting  switch  S2 
since  the  continuous  rating  is  exceeded.  This  cannot  be  achieved  electro- 
mechanlcally  since  the  diode  switch  can  fail  in  less  than  2  jxS  at  the  applied 
power  levels.  Probably  the  only  method  available  is  to  have  separate  drive 
circuitry  for  each  of  the  switches  in  Sj^. 
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5 . 3  Switching  Time 

The  "cherry  picker"  boom  system  (Ref.  1)  will  position  the  antenna 
at  a  maximum  height  of  40  feet  above  the  ground  for  a  vertical  incidence 
measurement.  The  signals  are  carried  from  the  equipment  truck  to  the  antenna 
by  coaxial  cable.  The  minimum  signal  return  time  is  thus  150  ns  referred  to 
the  duplexing  cir  ulator. 

The  timing  seqi  ;nce  is  shown  (Figure  5)  with  the  delay  times  indicated 
for  the  chosen  switches  .  If  the  switch  S2  were  mounted  on  the  boom  together 
with  the  transmitter  and  receiver  TWT's.the  signal  delay  between  switches 
would  become  important. 

5 . 4  The  Circulator 

The  circulator  must  be  able  to  pass  the  maximum  peak  power  of  20  watts 
and  have  sufficient  isolation  between  reverse  ports. 


6.0  FREQUENCY  CONVERSION 

In  order  to  reduce  the  bandwidth  of  the  receiver  it  is  necessary 
to  mix  the  received  signal  with  a  local  oscillator  signal  (fs  +  f»  )  displaced 
in  frequency  from  the  transmitted  signal  (fg)  by  the  chosen  IF  frequency 
(f^)  this  can  be  accomplished  in  three  ways. 

1.  A  double  frequency  conversion.  (Fig  .  6) 

2.  Using  a  voltage  tuned  filter.  (Fig.  7) 

3.  Using  quadrature  hybrid  phase  shifters  (Fig.  8) 

The  second  method  was  chosen  since  we  were  unable  to  obtain 
suitable  components  for  the  other  two  methods.  Mixers  were  not  available 
with  the  required  bandwidth  on  the  IF  port  for  method  1  and  Quadrature 
hybrid  phase  shifters  with  the  required  isolation  and  defined  phase  shift 
were  not  available  for  method  3. 

The  specifications  of  the  required  components  are  given  (Spec.  6)  . 
The  specifications  of  the  mixer  are  not  at  the  required  power  levels.  We 
must  therefore  estimate  the  conversion  loss  and  isolation. 
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We  will  examine  the  signal  levels  in  Lne  frequency  convert  ion 

system: 


Leveled  signal  from  PIN  modulator 

LO  power  required 

Conversion  loss 

Power  in  each  side  band 

YIG  insertion  loss 

Isolator  insertion  loss 

Effective  directional  coupler  insertion  loss 

Sideband  power  output 


=  15  mW 

=  100  mW  (estimate) 
=  -10  db  (estimate) 

=  1  . 5  mW 
=  5  db 
=  0.4  db 
=  1 .2  db 
=  0.33  mW 


The  IF  amplifier  mixer  combination  requires  a  local  oscillator 
power  of  greater  than  0.25  mW. 


7.0  LEVELING 

Two  closed  leveling  loops  are  included  in  the  system,  the  first 
to  level  the  transmitted  power  and  the  second  to  level  the  frequency  converter 
output.  The  method  of  leveling  is  shown  (Figure  9,  Ref.  4).  One  of  the 
required  amplifiers  is  included  in  the  sweep  frequency  oscillator;  the  other 
was  purchased  separately  (Reference  5) . 
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System  Power  Requirements 


Motor  Generator  Output 

15,000  watts 

Sweep  Oscillator 

350  watts 

TWT  Receiver 

8.5  watts 

TWT  Transmitter 

644  watts 

D.C.  Power  Supplies 

300  watts 

Oscilloscope  (585) 

550  watts 

Sampling  Oscilloscope 

500  watts 

Leveler  Amplifier 

2  watts 

Signal  Generator 

9  watts 

Tape  Punch 

138  watts 

Digital  Volt  Meter 

30  watts 

Antenna  Positioner 

1 , 300  watts 

Total  Load 

4,031  watts 
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Index  to  Figures 

1 .  Functional  arrangement  of  the  polypanchromatic  radar 

2.  Component  arrangement  in  system 

3.  Transmitter  TWT  with  associated  components 

4 .  Arrangements  of  duplexing  switches 

5.  Switching  timing  sequence 

6.  Methods  of  frequency  conversion 

7 .  Closed  loop  leveling 
a.  Duplexing  components 

9.  Power  rating  of  diode  switches 
10.  Frequency  spectrum  of  transmitted  pulses 


Antenna 


« 


Figure  1.  Slow- Sweep  System  Block  Diagram 


Van  Mounted  Boom  Mounted 

Equipment  Equipment 


* 


POLYPAN CHROMATIC  RADAR  -  RF  SECTION  BLOCK  DIAGRAM 


Transmitter 

T.W.T. 


Figure  4.  R.F.  Switch  Arrangement 


Switch  A2 


Switch  Aj 


Switch  A4 


Figure  5.  R.F.  Switching  Sequence 
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Fixed 

Frequency 

Source 


a  .  1st  IF  Below  Band 


b.  1st  IF  Above  Band 


Figure  6.  Frequency  Conversion  by  2-Stage  Mixing 
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Figure  7.  Frequency  Conversion  Using 
Voltage  Tuned  Filter 
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Figure  8  .  Frequency  Conversion  Using  Quadrature 
Hybrid  Phase  Shifter 


C22 


Unleveled 
Signal  In 


Leveled 
Signal  out 


Figure  9 


Circulator  Schematic 


Diode  Switch  3570  Schematic 


Diode  Switch  3603  Schematic 


Figure  10 


Incident  power  derating  factor 
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Module  case  temperature  — 
degrees  centigrade 


figure  11  .  Incident  Power  Rating  at  25°C  for  3603  at  4-8  GHz 
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Frequency  domain  amplitude  spectrum  of  a  pulse -modulated  carrier 


-*IT  |* 


* 


time 


Time  domain.  Envelope  of  pulse -modulated  carrier 


Overlapping  spectrum  due  to  slow  change  of  carrier  frequency 


Figure  12 
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Index  to  Specifications  of  Components 

1.  Transmitter  T  .W.T. 

2.  Sweep  Oscillator 

3 .  Receiver  T  .W  .T . 

4.  IF  Amplifier 

5.  Duplexing  Components 

6.  Frequency  Conversion  Components 
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SPECIFICATION  1 

Transmitter  T.W.T.  Amplifier 
Manufacturer  —  Alto  Scientific  Co.,  Inc. 
Type  —  20C  4.0  -  9.0  Q  60 
Frequency  —  4. 0-8.0  GHz 
Small  Signal  Gain  --  60  db 
Noise  Figure  --  30  db  nom.  35  db  max. 
Max.  R.F.  input  --  1  watt 
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Reference.  Alto  Scientific  Instruction  Manual  — Alto  Scientific  Engineering 
Note  5-66. 
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SPECIFICATION  2 

Sweep  Frequency  Oscillator 

Manufacturer  —  Hewlett  Packard 
Model  —  8690A  sweep  oscillator 
8693B  RF  unit 

Frequency  Range  — Adjustable  from  4-8  GHz  (manual  or  internal  sweep) 
Frequency  Accuracy  —  +  20  MHz 
Leveled  Power  Output  —  15  mW 

Sweeptime  Adjustable  in  4  decades  0.01  to  100  secs. 

Associated  Outputs  —  Sweep  Output 

Sweep  reference  output 

Blanking 

Market 

Leveling  Amp.  Input  —  -20  to  -350  mV 


Reference.  Hewlett  Packard  Manual  8690A,  8693B. 
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WATKINS  JOHNSON  COMPANY 
SPECIFICATION  3 

FINAL  DATA  SHEET 
LOW-NOISE  TWT  AMPLIFIER 


The  data  as  recorded  below  Is  the  result  of  tests  performed  by  our  test 
department  just  prior  to  shipment  of  this  tube,  and  Is  provided  as  a  special 
service  to  our  customers . 


Tube  Type  WJ-  423 


Serial  No.  5 


SUPERSEDES  ALL  PREVIOUS  DATA 


Frequency 

(GHz) 

VSWR  Maximum 

Small  Signal 
Gain  (dB) 

Noise  Figure 

(dB) 

* 

Saturation 
Power  Output 
(dBm) 

Input  Output 

4.0 

31.8 

7.3 

+  15.0 

4.4 

34.3 

6.8 

+  15.7 

4.8 

34.8 

6.7 

+  15.8 

5.2 

34.8 

6.5 

+  15.8 

5.6 

H  H 

34.6 

6.3 

+  15.7 

6.0 

•  •  •  • 

34.2 

6.0 

+  15.1 

6.4 

o  o 

33.5 

6.1 

+  14.7 

6.8 

•  • 

33.4 

6.2 

+  14.4 

7.2 

CM  CM 

32.4 

6.6 

+  13.4 

7.6 

A  A 

32.3 

6.9 

+  13.4 

8.0 

32.2 

7.3 

+  12.9 

Tube  Requirements 


*N.F.  measured  with  H.P.  noise  sources 
G347A  and  J347A  and  H.P.  noise  figure 
meter  model  340B. 


Input  Voltage:  105-125  V  ac 


Input  Current:  165  mA 

Input  Power:  8.5  Watts 


Tested  by:  M.  Slaughter _  Date:  21  March  1968 

3333  Hillvlew  Avenue  Stanford  Industrial  Park  Palo  Alto,  California 
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SPECIFICATION  4 
IF  Amplifier 

Manufacturer  —  R.H.G.  Labs 

Model  No.  Preamp.  MP  4-8m37,  Post  Amp.  EBT106MGC 
RF  Frequency  —  4-8  GHz 
IF  Frequency  —  250  MHz 
IF  Bandwidth  —  100  MHz 

Power  Gain  Preamp.  (RF  to  IF)  23  db,  Postamp.  55  db 
Noise  Power  —  4  db  (at  0.25  mW  LO) 

Maximum  Output  Preamp.  10  dbm,  Postamp.  0  dbm 


Reference.  RHG  Electrical  Test  Data  Form  MP1  and  Form  EVT1 
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SPECIFICATIONS  5 

Duplexing  Components  (Figure  10) 

Circulator  (Test  Data  Supplied) 

Manufacturer  —  Melabs  .  Model  6317-264 

Frequency  —  4-8  GHz 

Insertion  loss  port  1-2  —  0.36  db  max. 

1-3  —  19.0  db  max. 

Switch  Sj^  S4  (Two  each  in  series) 

Manufacturer  —  Hewlett  Packard 

Model  PIN  Diode  Switch  3570 

Frequency  —  4-8  GHz 

Insertion  loss  —  1 .5  db  max.  at  -20v  bias 

Isolation  —  38  db  min.  at  100  mA  bias 

VSWR  —  2.5:1  max.  at  -20v  bias 

Max.  RF  power  —  1  watt  at  25*  C. 

Maximum  bias  current  —  100  mA 

Maximum  bias  voltage - 36v 

Switching  speed  --  RF  ON  15  nS  max.  9  nS  typ. 

RF  OFF  10  nS  max.  6  nS  typ. 

Switching  speeds  are  measured  over  a  range  of  10  db  detected  power  (OFF) 
to  0.45  db  detected  RF  power  (ON) 

Switch  So 

Manufacturer  —  Hewlett  Packard 
Model  PIN  diode  switch  C00433603A 
Frequency  —  4-8  GHz 
Insertion  loss  —  1 .5  db  max.  Ov  bias 
Isolation  --  60  db  min.  -150  mA  bias 
VSWR  —  2.0:1  max.  at  Ov  bias 
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SPECIFICATION  5  CQNT'D. 

Max.  RF  power  —  see  Figure  11 

Max.  forward  bias - 675  mA 

Max.  reverse  bias - 100  v 

Switching  speed  —  RF  ON  50  nS  max. 

RF  OFF  40  nS  max. 
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SPECIFICATION  6 

Frequency  Conversion  Component 
Mixer 

Manufacturer  —  Sage 
Model  2543  NOCM 
Frequency  —  4-8  GHz 
IF  --  200  MHz 
Isolation  —  8  db  min. 

Conversion  los3  -  7.5  db  min.  (Sig.  -20  dbm  L.O.  3  dbm) 

Local  Oscillator 

Manufacturer  —  Frequency  Sources 
Frequency  --  250  MHz  (Mechanical  adj.) 

Power  —  200  mW  (Min.) 

Ylq  Filter 

Manufacturer  —  Watkins  Johnson 
Model  —  617-8  Four  Stage  Yig  Filter 

922-4  Voltage  to  Current  Driver,  250  MHz 
Frequency  —  4.25-8.25  GHz 
Bandwidth  —  45-80  GHz 

Tuning  Voltage  —  0-10  v.  • 

Insertion  Loss  —  5  db  max. 

Off  resonance  isolation  —  70  db  min. 

Off  resonance  spurious  —  50  db  min. 

Passband  ripple  —  1 .0  db 

Passband  ripple  spurious  —  1.5  db  (max.) 

Passband  VSWR  —  2:0  (max.) 

Limiting  level  —  10  dbm  min. 

Hysteresis  —  2-  MHz 

Deviation  from  linear  —  +8  MHz  max. 
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SPECIFICATION  6  CONT'D . 

Frequency  Deviation  —  JO  GHz  (0*  C  to  60*  C) 
Frequency  Selectivity  —  24  MHz/octave  (of  bandwidth) 
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A  DIGITAL  SYSTEM  TO  CONTROL  SAMPLING  MEASUF  MENTS 
FOR  A  POLYPANCHROMATIC  RADAR 


Its'hak  Dinstein 


Technical  Memorandum  133-1-3 


A  DIGITAL  SYSTEM  TO  CONTROL  SAMPLING  MEASUREMENTS 
FOR  A  POLYPANCHROMATIC  RADAR 

by 

Its'hak  Dlnstein 


The  digital  control  system  is  designed  to  control  the  time  sampling  in  the 
frequency  domain  of  reflected  pulses  in  a  wide  bandwidth  radar  system.  The 
memorandum  is  written  in  a  general  way  and  is  almost  independent  of  the  radar 
system.  The  block  diagram  of  the  whole  system  is  given  here  briefly. 


The  Digital  Control  Unit  will  be  operated  with  -5  Volts ,  (V  of  the  TTL's 

cc 

connected  to  ground  and  ground  of  the  TTL's  connected  to  -5  Volts)  in  order 
to  be  compatible  to  signals  to  and  from  the  sweep  OSC . 

The  format  of  this  memorandum  is  that  of  a  step-by-step  procedure. 
This  seems  to  describe  the  logical  and  switching  nature  of  this  system  better 
than  the  more  usual  narrative  description. 
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THE  PROBLEM 


1.1  Given  a  periodic  function  of  time  f(t)  (See  Figure  1) 

1.2  Divide  the  period  Tq  into  1000  increments. 

1.3  To  be  able  to  choose  any  time  interval  t^  to  t2  so  that  ^  ~  tl  =  Ti ' 
and  to  produce  an  output  signal  during  this  time. 

1.4  To  repeat  this  output  signal  during  the  time  t^  +  T  to  ^  +  T  (see 
Figure  1  )  and  again  during  the  time  t^  +  2t  to  t2  +  2t,  and  so  on 
n  times  till  the  time  t2  +  (n  -  1)  t. 

1 .5  Then  starting  again  from  tg,  to  choose  another  interval  t3  to  t4 

so  that  t4  -  t3  =  and  t3  -  tj  =  t2  (see  Figure  1 )  ,  and  to 

produce  an  output  signal  during  this  time  interval. 

1.6  To  repeat  this  output  signal  during  the  time  t3  +  -  to  t4  +  t  and 
so  on  n  times  till  the  time  t4  +  (n  -  1)t 

1.7  To  repeat  this  procedure  till  m  intervals  are  covered. 

THEORY  OF  OPERATION 

2.1  An  internal  clock  produces  a  time  base.  The  frequency  of  the  clock 

is  adjusted  to  be  (Figure  1 ) 

T0 

T0 

2.2  The  unit  of  time  is  equal  to  J^qq'  • 

2.3  The  timing  is  controlled  by  5  counters: 
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1)  Start  Counter  (STR.  CNT) 

2)  Bandwidth  Counter  (BND.  CNT) 

3)  Repetitions  Counter  (REP.  CNT) 

4)  Spacing  counter  (SPC  .  CNT) 

5)  Number  of  Bands  Counter  (NMB.  CNT) 

Each  counter  consists  of  3  decimal  digits.  The  counters  can 
be  preset  to  any  number  defined  by  the  DIGISWITCHES  (Figure  2) 

2.4  The  beginning  of  the  first  interval  of  time  is  defined  by  the  START 
COUNTER.  Tha  STR.  CNT  is  preset  to  the  complement  number  of 
tj  (Figure  1) .  Therefore,  the  counter  will  be  full  after  counting 
tj  pulses.  When  the  STR.  CNT  starts  to  count  at  t^,  it  reaches  tj 
when  it  is  full. 

2.5  The  length  of  the  time  interval  is  defined  by  the  BANDWIDTH  COUNTER. 
The  BND .  CNT  is  preset  to  the  complement  number  of  Tj .  (Figure 

1)  .  The  counter  will  be  full  after  counting  Tj  pulses  .  The  BND . 

CNT  starts  its  counting  when  the  STR.  CNT  starts  its  counting  when 
the  STR.  CNT  is  full,  therefore  it  will  be  full  at  t ^ . 

2.6  The  number  n  (see  sec.  1.4)  is  defined  by  the  REPETITION  COUNTER. 
The  REP.  CNT  is  preset  to  the  complement  of  n,  and  is  advanced  by 
one  after  every  interval  of  time.  When  the  REP  CNT  is  full,  the 
system  begins  to  cover  the  next  interval. 

2.7  The  beginning  of  the  next  interval  is  defined  by  the  SPACING 
COUNTER  and  by  the  START  COUNTER.  When  moving  to  the  Kth 
interval,  the  SPC.  CNT  will  count  (K  -  1)  x  pulses,  and  when  it 
will  be  full,  the  STR.  CNT  will  count  tj  pulses.  The  beginning 

of  the  Kth  interval  is  at  t  =  tj  +  (K  -  1)  x  r2 . 
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2.8  The  sequency  is  finished  when  the  NUMBER  OF  BANDS  COUNTER 

is  full.  The  NMB.  CNT  is  preset  to  the  complement  number  of 
m  \see  Sec.  1.7).  The  counter  is  advanced  by  one  when  the 
system  is  moving  from  one  interval  to  the  next. 

OPERATION 

3 . 1  Calibration 

3.1 .1  The  CALIBRATION  COUNTER  (CLB.  CNT)  starts  counting  at  tQ, 
triggered  by  an  external  start  marker. 

3.1.2  When  the  CLB-  CNT  counts  999  pulses  this  state  is  decoded  and 

then  the  D  input  of  FF  10  (see  the  overall  scheme)  is  going  high. 

If  the  external  end  marker  of  Tn  (Figure  1  )  comes  between  the 
th  " 

999  pulse  and  the  1000  pulse,  the  clock  is  calibrated  and  the 
calibration  lamp  will  be  on. 

3.1.3  If  the  CLB.  CNT  is  full  before  Tq,  FF  11  will  be  triggered  and  the 
overflow  lamp  will  be  on.  It  indicates  that  the  frequency  of  the 
clock  should  be  deer  :ased. 

3.1.4  When  the  clock  is  calibrated,  the  CAL  lamp  is  on  and  the  overflow 
lamp  is  off. 

3.2  The  RUN  MODE  (see  overall  scheme) 

3.2.1  All  the  Digiswitches  have  to  be  adjusted  to  the  proper  numbers. 

The  system  has  to  be  reset  and  then  preset. 

3.2.2  The  given  function  (Figure  1  is  produced  by  the  sweep  OSC.  The 
sweep  OSC  has  to  be  triggered  manually,  to  start  the  sequence. 
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3.2.3  When  the  sweep  begins,  FF  1  is  triggered  to  on  and  then: 

a)  The  STR.  CNT  begins  to  count  cr-pulses  (through  gate  6) 

b)  FF  3  is  reset. 

3.2.4  When  the  STR.  CNT  is  full,  it  triggers  FF  7  and  gate  19  is  opened, 
then: 

a)  The  SAMPLER  and  the  RF  UNIT  DRIVERS  are  triggered  by  a 
pulses  (this  is  the  output  signal). 

b)  The  BND.  CNT  is  not  counting  a  pulses  (through  gate  19) 

3.2.5  When  the  BND.  CNT  is  full,  it  triggers  FF  8,  and  gate  19  is 
closed . 

3.2.6  When  the  sweep  comes  to  its  end: 

a)  FF  7  and  FF  8  are  reset  by  p  pulses  (through  gate  3  and  gate  8) . 

b)  The  STR.  CNT  and  BND.  CNT  are  preset  by  the  p-pulses. 

c)  The  REP.  CNT  is  triggered  (through  gate  3)  and  it  is  advanced 
by  one . 

3.2.7  If  the  REP.  CNT  is  not  full,  the  sweep  OSC  is  triggering  itself 
(by  the  blanking  pulse  through  gate  1)  and  the  procedure  repeats 
again. 

3.2.8  When  the  REP.  CNT  is  full,  then: 

a)  It  triggers  the  monostable,  and  gate  1  is  closed  and  the  sweep 
OSC  is  not  triggered  by  itself. 

b)  An  output  signal  from  the  monostable  indicates  that  the  system 
is  moving  to  the  next  Interval. 

c)  FF  2  is  triggered  by  the  REP.  CNT,  gate  7  is  opened,  the  number 
of  SPACES  CNT  preset  by  p-pulses  and  the  NMB .  CNT  is  advanced 

by  one  by  p  pulse  (through  gate  7) . 

3.2.9  When  the  monostable  returns  to  its  stable  state  it  triggers  the  sweep 
OSC. 
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3.2.10  The  comparator  is  now  not  equal,  because  the  number  of  SPACES 
CNT  was  preset  and  the  NMB.  CNT  was  advanced  by  one.  There¬ 
fore  ,  gate  5  is  now  open  and  the  SPC .  CNT  is  counting  a  pulses . 

3.2.11  When  the  SPC .  CNT  is  full: 

a)  It  advances  the  number  of  SPACES  CNT  by  one 

b)  It  triggers  FF  4,  then  the  p  pulses  preset  the  SPC .  CNT  and  the 
next  a  pulses  reset  FF  4. 

3.2.12  The  comparator  is  now  equal  and  therefore  gate  5  is  closed  and 
gate  6  is  opened. 

3.2.13  The  STR.  CNT  counts  a-pulses  and  the  procedure  continues  as 
before . 

3.2.14  When  the  NMB.  CNT  is  full,  the  sequence  is  completed,  FF  9 
is  on,  gate  19  is  closed  and  there  is  no  more  output  signal. 


DECADE  CO 


LOGIC  CONTROL  SYSTEM  FUNCTIONAL  DIAGRAM 


Appendix  E 


MULTIPLE  IMAGE  CLASSIFICATION  WITH  SPATIAL 
AND  MEASUREMENT  SPACE  CLUSTERING 


R.  M.  Haralick 


Technical  Memorandum  133-6 


MULTIPLE  IMAGE  CLASSIFICATION  WITH  SPATIAL 
AND  MEASUREMENT  SPACE  CLUSTERING 


R.  M.  Haralick 
Remote  Sensing  Laboratory 
Center  for  Research  in  Engineering  Science 
University  of  Kansas 
Lawrence,  Kansas  66044 


We  have  done  work  in  texture  analysis,  adaptive  •  ustering,  and 
spatial  clustering.  We  have  also  implemented  various  optimizing,  sta¬ 
tistical  and  Bayes  decision  theory  techniques  for  the  computer.  Included 
in  this  report  is  a  preliminary  description  of  this  research.  Detailed  tech¬ 
nical  reports  will  be  forthcoming  for  each  area .  The  spatial  clustering 
report  will  be  out  shortly  and  it  is  not  discussed  here. 

Basic  research  in  texture  analysis  (as  part  of  the  color  combination 
and  texture  study  of  multiple  images)  was  undertaken  because  we  feel 
that  the  local  spatial  structure  of  the  image  densities  has  much  informa¬ 
tion  about  its  classification.  This  has  been  further  illustrated  with  our 
spatial  clustering  results.  On  a  test  set  of  three  images,  the  classifica¬ 
tion  obtained  using  spatial  clustering  had  about  half  the  errors  which  the 
classification  using  measurement  space  clustering  had. 

The  adaptive  clustering  technique  which  was  briefly  described 
in  the  study  plan  was  programmed.  Difficulty  was  found  in  getting  the 
procedure  to  converge  for  arbitrary  starting  points.  If  the  starting  point 
was  sufficiently  close  to  a  solution  then  the  procedure  worked  fine. 
Attempts  are  now  being  made  to  simplify  the  algorithm  by  breaking  it  up 
into  parts  which  work  independently.  In  this  manner  each  part  can  con¬ 
verge  without  upsetting  the  other  parts.  We  are  testing  various  measures 
of  similarity  to  find  the  one  which  works  best. 

We  are  also  investigating  the  possibility  of  combining  the  adap¬ 
tive  techniques  with  standard  optimizing  procedures.  In  this  vein  we 
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have  written  gradient  search,  dynamic  allocation,  generalized  matrix 
inverse,  matrix  inverse  accuracy  increaser,  and  eignvector-eigenvalue 
routines.  We  have  continued  our  work  with  Bayes  theory,  updating  our 
Bayes  program  to  include  many  general  options.  A  brief  description  of 
both  Bayes  theory  and  the  programs  which  can  be  used  in  optimizing  pro¬ 
cedures  are  included. 


COLOR  COMBINATION  AND  TEXTURE  STUDY  OF  MULTIPLE  IMAGES 

The  intent  of  this  task  was  to  examine  color  and  texture  combina¬ 
tions  of  imagery  during  the  summer  of  1968.  However,  it  became  apparent, 
in  our  discussions  with  the  visual  perception  psychologist  Dr.  Epstein, 
that  if  we  wanted  to  study  both  color  and  texture  together,  we  would  first 
have  to  study  texture  in  black  and  white.  Much  more  is  known  about  color 
than  about  texture,  and  since  we  eventually  want  to  perform  a  factorial 
experiment  involving  both  color  and  texture  we  must  be  equally  knowledge¬ 
able  about  both.  Our  first  goal  was  to  create  a  set  of  abstract  textures 
with  different  texture  gradients  in  order  to  determine  where  subjects  would 
perceive  the  most  convincing  edge  or  boundary.  With  this  information  we 
hoped  to  be  able  to  create,  by  superposition,  multiple  texture  images 
which  were  simple  textural-map  representations  of  actual  radar  imagery. 

Our  first  problem  was  to  decide  whether  dots,  lines,  or  other  vari¬ 
ous  geometric  figures  should  be  the  basic  pattern  in  the  texture.  We 
decided  upon  short  lines  because  groups  of  scattered  short  lines  had 
least  similarity  to  common  textured  objects:  other  texture  patterns  such  as 
various  types  of  checkerboard  patterns  could  easily  be  perceived  as  a 
wall,  floor  or  receeding  surface;  it  was  this  kind  of  identification  which 
we  wished  to  avoid.  In  addition,  textured  patterns  of  lines  are  easy  to 
generate  on  the  plotter  available  at  the  KU  Computation  Center. 

The  principal  goal  was  to  produce  a  continuum  of  computer  gen¬ 
erated  textures  which,  as  a  function  of  known  parameters,  range  from 
those  suggesting  a  uniform  flat  surface  to  those  strongly  suggesting  two 
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surfaces  meeting  at  a  boundary,  thus  forming  an  edge.  Two  different 
approaches  have  been  tried:  one  using  a  simple  geometric  projection  model 
to  determine  the  distribution  gradients,  and  the  other  using  a  set  of  gradients, 
resulting  from  no  single  model,  but  empirically  showing  the  desired  visual 
effect. 

Boundaries  are  caused  by  an  abrupt  change  of  texture  gradients. 

An  edge  is  produced  by  introducing  a  discontinuity  of  the  gradient  of  the 
line  length  and  line  position  distributions.  Visual  displays  were  made  by 
Joining  two  gradients  in  one  texture  pattern.  At  the  edge  (where  the  grad¬ 
ients  interfaced)  the  distributions  were  continuous,  so  that  only  the  change 
in  gradients  could  be  observed.  The  two  texture  patterns  used  were  as 
follows: 

The  first  texture  pattern  was  produced  by  simulating  the  image  on 
the  retina  occuring  when  a  person  sees  a  receeding  surface  met  by  a  wall. 

One  gradient  was  the  null  gradient,  that  is,  the  distributions  of  length  and 
position  were  uniform.  The  second  gradient  was  that  of  a  non- vertical 
plane  surface.  To  determine  the  resultant  distribution  on  the  retina  image 
it  was  necessary  only  to  compute  the  projection  of  a  uniform  distribution 
on  a  plane  onto  a  spherical  surface . 

The  second  set  of  texture  patterns  was  constructed  arbitrarily  using 
''reasonable"  analytic  gradients.  The  functions  are  of  the  form 

f(y)  *  c(i-b)m/2 
where  c  is  a  constant 

b  is  a  uniformly  distributed  random  variable  over  (-.5,  +  .5) 

An  extensive  search  of  the  psycho- physical  literature  was  under¬ 
taken  to  determine  what  was  the  best  experimental  methodology.  Under 
consideration  were  the  following  methods:  paired  comparison,  multiple 
comparison,  rank  order  and  absolute  value.  In  the  paired  comparison 
methods,  two  texture  gradients  at  a  time  are  presented  to  each  subject, 
and  the  subject  must  choose  which  one  has  the  most  convincing  edge. 

In  the  multiple  comparison  method  three  or  more  texture  gradients  are 
presented  at  one  time,  and  the  subject  must  choose  which  one  has  the 
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most  convincing  edge.  In  the  rank  order  method,  all  the  texture  gradients 
are  presented  at  one  time  and  the  subject  must  order  the  gradients  from  the 
one  with  the  most  convincing  edge  to  the  least  convincing  edge.  In  the 
absolute  value  method,  the  subject  is  presented  with  one  texture  gradient 
at  a  time,  and  he  must  assign  some  integer  between  0  and  12  for  each  one; 

12  indicates  no  edge  and  0  indicates  an  absolutely  sharp  edge. 

The  absolute  value  method  was  chosen  because  it  allows  an  aver¬ 
age  numerical  quantity  to  be  calculated  for  each  texture  gradient;  from 
the  numerical  quantity  the  gradients  cannot  only  be  rank  ordered  but  the 
distance  from  one  gradient  to  the  next  may  be  estimated.  The  other  meth¬ 
ods  assume  nothing  and  only  allow  a  rank  ordering  to  be  made.  The  dis¬ 
advantage  of  the  absolute  value  method  is  that  it  assumes  the  subject  will 
unconsciously  use  a  ratio  scale  to  assign  numbers  (the  number  assigned 
to  a  texture  gradient  which  has  an  edge  twice  as  convincing  as  another, 
will  be  assigned  a  number  which  is  twice  as  great  as  the  other). 

The  computer  programs  were  written  to  generate  the  texture  gradients. 
Difficulty,  however,  was  experienced  in  plotting  the  lines.  The  plotter 
was  broken  for  days  and  weeks  at  a  time;  it  was  working  less  than  fifty 
per  cent  of  the  time.  By  the  middle  of  August  a  satisfactory  model  for 
generation  of  the  texture  gradients  was  found  and  plots  had  been  made. 

There  was  enough  time  to  set  up  and  do  the  first  experiment  before 
Dr.  Epstein  left  the  University. 

Presently,  line  texture  displays  are  being  prepared  with  various 
line  orientations.  Plans  include  the  production  of  textures  which  have 
basic  elements  which  enclose  larger  areas  than  short  lines  do.  We  are 
also  investigating  the  possibility  of  using  a  line  printer  to  generate  the 
textural  displays. 


PROGRAMS  USED  IN  OPTIMIZING 


Matrix  Inversion 

1 .  The  first  of  the  programs  dealing  with  matrix  inversion  is  a 
subroutine  called  GINV.  This  program  calculates  the  generalized,  or 
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pseudo- inverse  of  a  rectangular  matrix  A  which  may  be  singular  or  non¬ 
singular.  The  psuedo  inverse  of  A  can  be  defined  as  follows:  Given  the 
matrix  A  of  order  N  x  M,  the  generalized  inverse  is  the  matrix  A  of  order 
M  x  N  which  satisfies  the  equation 

AA'A  =  A.  (1) 

2.  The  second  of  these  routines  dealing  with  matrix  inversion  is 
used  to  increase  the  accuracy  of  the  inverse  of  a  non- singular,  square  ma¬ 
trix.  Given  the  matrix  A  and  its  estimated  inverse  B,  subroutine  ACCINC 
calculates  a  new  inverse  B*  whose  accuracy  is  increased  by  an  arbitrary 
degree  specified  by  an  input  parameter,  EPS.  The  success  of  ACCINC 
requires  that  one  condition  be  fulfilled.  This  condition  is  described  in 
the  equation  below. 

|  1 1  -  AB  |  |  <1  (2) 

where  I  is  the  identity  matrix. 

The  calculation  process  of  GINV  requires  that  the  inverse  of  a  non¬ 
singular  square  matrix  be  found,  and  consequently  the  user  may  desire  a 
higher  degree  of  accuracy  than  is  normally  provided  by  many  standard 
Inversion  routines.  Therefore  subroutine  ACCINC  can  be  incorporated  into 
GINV  by  exercising  one  of  the  user's  options  of  GINV. 

Eigenvalues  and  Eigenvectors 

1 .  Subroutine  BISECT  was  written  to  calculate  the  eigenvalues 

of  a  symmetric  tridiagonal  matrix  by  the  method  of  bisection.  The  original 
program  was  written  by  W.  Borth,  R.  S.  Martin,  and  J.  H.  Wilkinson  in 
ALGOL  and  represents  one  of  the  fastest  and  most  recent  variations  of  the 
bisection  method.  It  was  converted  to  FORTRAN  IV  and  put  in  subroutine 
form.  Later  it  was  combined  with  the  following  two  routines  to  form  a 
package  to  calculate  eigenvalues  and  eigenvectors. 

2.  Subroutine  TRIDI  was  derived  after  testing  many  existing  rou¬ 
tines  which  found  the  tridiagonalized  form  of  a  symmetric  matrix.  TRIDI 
embodies  the  desirable  characteristics  of  all  of  these  routines  including 
efficiency  in  calculation,  increased  accuracy  and  minimum  storage  require- 


E5 


merits.  As  mentioned  before,  this  was  included  in  the  eigenvalue- eigen¬ 
vector  package. 

3.  The  last  of  these  programs  was  developed  in  much  the  same 
way  as  TRIDI.  Subroutine  VECTOR  is  used  to  calculate  the  eigenvectors 
of  the  matrix  and  incorporates  the  best  characteristics  of  the  many  routines 
that  were  tested. 

Optimal  Search  Techniques 

1 .  The  first  of  the  programs  dealing  with  optimal  search  techniques 
applies  itself  to  the  one-dimensional  allocation  problem.  The  routine  uses 
dynamic  programming  to  maximize  the  non-linear  function, 

R(X1,X2,...,Xn)  =  g^Xj)  +  g2(X2)  +  ...+gn(Xn)  (3) 

where  g^  is  the  i^  activity  function  and  X^  is  the  allocation  to  that  act¬ 
ivity.  This  is  solved  over  all  those  points  in  a  region 

^>0  1=1,2 . N  (4) 

which  also  satisfy  the  constraint 
N 

S  X.  <  XMAX  (5) 

i=l  1 

where  XMAX  is  the  total  resource. 

Since  the  solution  of  this  is  done  in  an  incremental  fashion,  the  g^  's  need 
not  be  continuously  defined  between  0  and  XMAX.  Rather  they  need  only 
to  be  defined  for  a  finite  set  of  points  each  separated  from  the  next  by  a 
delta  defined  by 

DELTA  =  XMAX/NINT  (6) 

where  NINT  is  the  number  of  intervals  into  which  the  range  between  0  and 
XMAX  is  divided  into.  Both  XMAX  and  NINT  are  input  variables  to  subrou¬ 
tine  DYNMC. 

2.  If  we  are  interested  in  maximizing  the  sum  of  a  set  of  N  acti¬ 
vities  which  are  functions  of  M-variables  (M  =  2,  3, . . .),  we  wish  to 
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solve  an  M-dimensional  allocation  problem.  Similar  to  the  previous  pro¬ 
blem,  we  wish  to  maximize  the  non-linear  function, 


KfV-V . V;X12-X22' 

2  g.  (X.1,  X.2,...,X.M) 
i=l  1  1  1  1 


•*nM> 

(7) 


where  the  superscript  indicates  the  dimension  number  in  an  M-dimensional 
space.  This  is  subject  to  the  conditions 


2  X. j  <  XMAXJ  X.J  >  0  i  =  1 , . . . ,  N;  j  =  1 , . . . ,  M 
i=l  1  1 

For  a  two  dimensional  problem,  (7)  becomes  the  following: 


R(X 


1'  2 ' 


,XN' 


P  12‘ 


-yn>- 


N 

2 

i=l 


(8) 


(9) 


in  which  the  superscript  notation  has  been  dropped.  By  keeping  the  X 
allocation  vector  constant,  (9)  reduces  to  the  same  one- dimensional 
allocation  problem  as  described  previously. 

N 

R(V  Y2 . V  -  =  «i<VYi>  <10> 

Therefore  if  we  supply  an  initial  guess  of  the  X  allocation  vector,  an 
approximate  solution  of  the  Y  allocation  vector  can  be  obtained  by  using 
the  first  program  described  In  this  section.  Further  by  using  the  Y  solu¬ 
tion  as  the  constant  vector,  an  approximate  solution  of  the  X  allocation 
vector  can  now  be  found.  Repeating  this  procedure  until  the  solutions 
remain  unchanged  for  two  successive  approximations,  a  maximum  of  R  is 
obtained.  However,  this  may  not  be  the  absolute  maximum  of  R  but  only 
a  relative  maximum.  Finding  one  of  the  absolute  maximums  may  be  highly 
dependent  on  the  initial  guess  because  of  the  nature  of  the  activities. 
Therefore,  more  than  one  attempt  should  be  made  in  finding  the  absolute 
maximum  each  time  using  a  different  initial  guess. 
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To  apply  this  method  of  successive  approximations  to  the  more 
general  case  of  m- dimensions,  m-1  dimensions  need  to  remain  constant 
to  reduce  it  to  the  one-dimensional  allocation  problem.  Therefore,  m-1 
initial  guess  vectors  should  be  suppled  to  start  the  calculation  process. 

The  two  subroutines  that  were  developed  to  handle  this  problem 
were  called  SUCESN  and  DYNMCN  where  SUCESN  controls  the  successive 
approximation  methods  and  DYNMCN  solves  the  one- dimensional  alloca¬ 
tion  problem. 

3.  Subroutine  STEEP  was  written  to  solve  for  the  minimum  of  a 
unimodal  function  of  N  variables  in  a  bounded  region  by  the  method  of 
steepest  descent.  If  the  function  is  represented  by 

g(Xj,  x2,...,xN)  (11) 

subject  to  the  conditions 

di  <  xi  <  bi  i  =  1 »  2, . . .  ,N  (12) 

then  subroutine  STEEP  will  solve  for  the  true  minimum  of  g  after  being 
supplied  with  an  initial  guess  of  this  minimum. 

The  first  step  of  the  solution  is  to  find  the  gradient  vector  at  the 
initial  point.  Next  a  search  is  made  in  a  direction  opposite  to  the  gradient 
vector  to  find  a  minimum  between  the  initial  point  and  the  nearest  boundary. 
This  is  accomplished  by  using  a  Fibonacci  search  procedure.  When  the 
minimum  is  found,  then  this  can  be  used  as  the  new  initial  point  and  the 
above  steps  can  be  repeated  to  find  a  new  minimum.  This  continues  until 
the  minimum  remains  unchanged  for  two  successive  searches  which  means 
that  the  true  minimum  has  been  found. 
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THE  BAYESIAN  APPROACH  TO  IDENTIFICATION 
OF  A  REMOTELY  SENSED  ENVIRONMENT 


Using  remote  sensors  we  can  make  measurements  of  an  environment. 
The  set  of  measurements  made  will  be  called  the  data  set.  Our  job  is  to 
examine  the  data  set  in  order  to  identify  what  the  environment  is  made  up  of: 
our  problem  is  how  should  we  do  it?  In  what  follows  we  describe  the 
Bayesian  decision  approach  with  a  deterministic  decision  rule. 

We  assume  that  distinct  boundaries  enclose  a  limited  environment, 
which  is  made  up  of  small-area  patches,  one  next  to  the  other.  The  identi¬ 
fication  of  the  environment  consists  of  identifying  each  small-area  patch 
within  one  category  of  a  given  set  of  categories.  We  assume  that  such  an 
identification  is  sensible  and  possible. 

In  order  to  make  any  identification  we  must  have  knowledge  con¬ 
cerning  which  kind  of  measurements  are  typical  measurements  of  the  cate¬ 
gories  we  wish  to  identify.  This  knowledge  is  succinctly  contained  in  a 
classification,  which  is  a  mapping,  associating  with  each  measurement 
the  category  to  which  it  is  most  typical  —  given  a  specific  decision  cri¬ 
terion.  Therefore,  if  we  are  to  identify  measurements  in  a  data  set  we 
must  have  a  classification. 

How  do  we  obtain  a  classification?  We  perform  an  information 
gathering  experiment.  From  the  population  of  all  environments,  we  sample 
one  or  a  few  in  which  it  is  possible  to  identify  many  small-area  patches  within 
each  category  of  interest.  The  proportion  of  occurrence  of  each  category  in 
the  sampled  environment (s)  does  not  have  to  be  representative  of  the  average 
probability  of  occurrence  of  each  category  in  the  entire  environmental  popu¬ 
lation.  However,  if  we  have  no  information  regarding  the  average  probability 
of  occurrence  of  each  categi  y  in  the  environmental  population,  then  we 
would  want  to  choose  the  sampled  environment(s)  so  that  the  proportion  of 
occurrence  of  each  category  in  the  sampled  environment(s)  is  an  unbiased 
estimate  of  the  average  probability  of  occurrence  of  each  category  in  the 
environmental  population.  In  either  case,  the  small-area  patches  within 
each  of  the  sampled  environment(s)  do  have  to  be  representative  of  the  cate¬ 
gories  with  which  they  are  identified. 
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With  each  of  our  sensors,  we  measure  each  small-area  patch  in  the 
chosen  environments.  From  photo- interpretation  or  field  studies,  the  en¬ 
vironments  are  examined  first  hand,  and  an  identification  of  each  small- 
area  patch  is  made.  The  sequence  of  such  identifications  is  called  the 
"ground  truth  identification"  or  simply  "ground  truth."  It  is  from  the  data 
set  (the  sequence  of  measurements)  and  the  ground  truth  (the  sequence  of 
identifications)  that  we  can  find  a  Bayes  classification. 

At  this  point  we  must  introduce  some  mathematical  notation. 


Let  C  =  { c ^  be  the  set  of  K  given  categories;  c.  is  the  symbol  used  for 
i=l  1 

th 

the  i  category.  We  suppose,  for  convenience,  that  each  sensor  produces 

only  one  number  for  each  measurement  it  makes  of  a  small-area  patch.  We 

,th 


suppose  further,  that  the  j  sensor  must  produce  a  number  belonging  to  its 
range  set  =  { 1^  ,  1^,  •  •  • ,  }  •  This  supposition  is  fully  in  accord 

with  reality,  since  the  output  of  any  sensor  is  always  equivalent  to  a  pointer¬ 
reading  on  a  dial.  Pointer-readings  can  never  be  discerned  precisely,  and 

th 

are  thus  discerned  approximately  to  third,  or  fourth,  or,  . .  .  ,N  place 
accuracy. 

Measurement  space  M  is  the  set  of  all  measurements  which  are 
possible  to  make  with  the  set  of  S  sensors.  M  is  conveniently  described 
as  the  cartesan  product  of  the  range  sets;  M  =  Lj  x  L2  x  ...  x  Lg.  This 
is  the  set  of  measurements  which  contain  for  elements,  all  the  possible 
numbers  produced  by  sensor  one,  combined  with  all  the  possible  numbers 
produced  by  sensor  two,  . . . ,  combined  with  all  the  possible  numbers  pro¬ 
duced  by  sensor  S.  For  convenience  we  number  the  measurements  in  M; 


M = {m  }N , 
n=l 


where  N  is  the  total  number  of  elements  in  measurement  space. 


Finally  we  must  provide  a  goodness  criterion;  thus,  we  introduce  a  gain 

function  g.  g(c.,c.)  is  our  economic  gain  if  we  identify  a  measurement  as 
1  J  th 

belonging  within  the  i  category  when  that  measurement  was  made  of  a 

th 

small-area  patch  actually  belonging  within  the  j  n  category. 
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We  have  already  mentioned  that  a  classification  is  a  mapping  or 
rule  which  associates  with  each  measurement  mn  in  M,  the  category  c^, 
to  which  it  is  most  typical  —  according  to  some  decision  criterion.  Our 
decision  criterion  is  economic;  "most  typical  to"  translates  to,  "that 
association  by  which  we,  on  the  average,  gain  the  most  economically." 
Therefore,  according  to  our  decision  criterion,  we  can  judge  each  possible 
classification.  That  classification  which  enables  us  to  gain  the  most,  on 
the  average,  is  the  classification  which  is  best;  it  is  that  classification 
which  we  wish  to  find. 

Let  us  now  examine  how  the  average  gain  may  be  calculated.  Let 
f  be  a  classification  mapping,  f  is  a  function  whose  domain  is  the  set 
M,  and  whose  range  is  the  set  C;  f;  M  —  C.  For  each  element  mn  e  M 
the  function  associates  one  and  only  one  category  ci  e  C.  We  define  the 
characteristic  function  h^  for  f  as  follows:  for  every  mn  e  M,  c^  e  C, 

hf(Cj,m  )  =  1  if  and  only  if  f(mn)  -  ci 
0  otherwise 

In  other  words  h^c^m^)  is  1  if  and  only  if  the  classification  f  identifies 
the  measurement  mn  as  belonging  within  the  category  c. .  The  average 
gain  A  for  the  classification  f  is  easily  seen  to  be: 

K  K  N 

M£,=  III  9(ck'cl)hf(0lc'mn)l,(nin|c1)P(cl) 

i=l  k=l  n=l 


where  P(mn | c^)  is  the  conditional  probability  that  the  measurement  m  will 
be  made  of  a  small-area  patch  given  that  the  patch  belongs  within  category 
Cj#  PCc^  Is  the  probability  that  any  small-area  patch  of  the  environments 
In  the  population  belongs  within  category  c^,  and  g(c^ ,  c^)  is  the  amount 
gained  If  a  patch  which  actually  belongs  within  category  c^  is  identified 
within  category  c^. 
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Of  the  four  terms  in  the  summation,  g(c^,cp  is  specified  as  part  of 
the  identification  goodness  criteria,  h^(c^,mn)  is  defined  from  the  classifi¬ 
cation  f,  P(mn|Ci)  will  be  determined  from  the  data  gathered  in  the  experi¬ 
ment,  and  P(cJ  is  an  additional  _a  priori  probability  which  we  will  have  to 
specify.  Let  us  now  examine  in  detail  how  the  conditional  probabilities 
are  determined  from  the  experimental  data. 

The  data  set  is  a  sequence  D  of  R  measurements; 

D  =  <(  m  ,m  ,...,m  \  . 

rl  r2  rR ' 

The  ground  truth  corresponding  to  sequence  D  is  a  sequence  T  of  R  not 

necessarily  different  category  identifications;  T  =  <fc  ,c  ,...,c  \  . 

rl  r2  rR 

Let  #  be  the  counting  measure.  #(D)  is  the  number  of  elements  in  the 
sequence  D;  thus,  #(D)  =  R.  A  sequence  is  really  a  function  whose  domain 
is  the  set  of  integers  I.  The  data  set  D  is  then  a  function  which  associates 
with  each  integer,  a  measurement;  D:  I  -*  M.  The  ground  truth  T  is  also  a 
function  and  it  associates  with  each  integer  a  category;  T:  I  -*  C.  D(7),  for 
example,  is  then  just  the  seventh  element  in  the  sequence  D;  D(7)  =  m  . 

_!  r7 

D  (m)  is  the  set  of  all  integers  i  for  which  D(i)  =  m.  The  statistic 

P (mn I c i)  estimating  P(mn|c.)  is  defined  as 

a  ,  ®'1(mJAT'1(cil)  . 

P(m  c.)  =  - » -  when  #  (T  1  (c.))  /  0 

n  1  #(T'1(ci))  1 

=  0  otherwise 

P (mn | c^)  is  the  number  of  integers  which  are  associated  with  the  measure¬ 
ment  mn  in  the  sequence  D  and  with  the  category  c^  in  the  sequence  T, 
divided  by  the  number  of  integers  associated  with  the  category  Cj  in 
sequence  T.  Stated  simply,  P(mn|c.)  is  just  the  number  of  times  the  mea¬ 
surement  mn  was  made  of  a  small-area  patch  belonging  within  category  c^, 
divided  by  the  number  of  times  a  small-area  patch  belonged  within  the  cate¬ 


gory  Cp 
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The  a  priori  probabilities  P(cp  can  either  be  estimated  from  the 
sampled  data  set  (if  this  data  set  is  representative  of  the  population)  or 
from  our  foreknowledge  of  the  population  of  environments.  If  we  can  assume 
that  the  few  environments  we  have  chosen  to  sample  for  our  experiment  are 
representative  of  the  population,  then 


P(cA)  = 


#(T_1  (c^) 
R 


is  a  reaonable  estimate.  If  we  cannot  make  such  an  assumption  and  we 
believe  that  a  small-area  patch  is  just  as  likely  to  belong  within  one  cate¬ 
gory  as  within  another,  then  P(c.)  =  1/K  is  a  reasonable  estimate. 

A  1  A 

From  the  estimates  P(m  c.)  and  P(c.)  we  may  estimate  the  average 
a  n 1  i  i 

gain  A  for  any  classification  f.  As  before  let  h^  be  the  characteristic 
function  for  f. 


hf(ck,mn)  =  1  if  and  only  if  f(mn)  =  ^ 
=  0  otherwise 


K  K  N 

X(f)=I  I  I  fltekiel)hf(c|t.mn)P(mn|c1)f(c1). 
1=1  k=l  n=l 


A 

We  seek  the  Bayes  classification  f*  which  maximizes  A.  f*  is 

easily  defined.  For  each  measurement  mn  and  for  any  classification  f, 

there  will  be  one  and  only  one  category  c.  such  that  h,(c.,m  )  =  1.  Con- 
/\  I  *  J  ^ 

sider  the  amount  a(Cj,mn)  gained  due  to  the  identification  of  measurement 

m  as  belonging  within  category  c,. 
n  j 

K 

a(Cj,mn)  =  ^  gte^c^PfmJc^PfCj) 
i=l 


I 
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The  maximum  A(f)  is  certainly  achieved  if  for  each  measurement  m  , 

/V  *  * 

f(mn)  =  Cj  where  c^  maximizes  a(Cj,mn).  Therefore  we  just  have  to  compute 
a(Cj,mn)  for  j  =  1 ,  2, ...  ,K  to  determine  which  category,  Cj,  maximizes  it. 
Then  we  define  f*(mn)  =  c^. 

In  this  manner  we  can  define  how  to  best  identify  each  measurement 
which  actually  occurred  in  the  data  sequence  D.  However,  there  may  be 
many  measurements  in  measurement  space  M  which  did  not  occur  in  the 
data  sequence.  How  should  these  measurements  be  identified  in  the 
classification?  Since  we  have  no  data  or  statistics  for  these  measurements 
it  seems  that  we  have  no  way  to  deal  with  them!  Here  we  must  draw  upon 
our  knowledge  of  the  structure  of  reality.  We  know  that  in  any  environment, 
if  a  measurement  m  is  made  of  a  small-area  patch  belonging  within  category 
c^  then  it  is  likely  to  make  measurements  m  +  6  for  other  small-area 
patches  which  also  belong  within  category  c^.  If  a  measurement  m  is 
typical  of  category  c^,  then  for  small  6,  m  +  6  is  also  typical  of  category 
Cj.  Similar  or  close  measurements  are  usually  associated  with  similar  or 
the  same  categories.  Thus  in  the  classification  we  can  identify  a  measure¬ 
ment  m,  which  did  not  occur  in  the  data  sequence,  with  the  category 
associated  with  m' ,  its  nearest  neighbor. 

The  part  of  the  classification  f*  which  was  defined  by  means  of 
the  statistics  generated  by  the  experiment  is  called  a  Bayes  Classification 
and  hence  the  name  "Bayesian  approach."  The  part  of  the  classification 
which  is  not  Bayesian  is  said  to  be  defined  by  a  nearest  neighbor  search. 
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IDECS*  SYSTEM  DEVELOPMENT 
G.  Kelly,  G.  Dalke 


Introduction 


The  IDECS  System  consists  of  real-time  and  near  real-time  image 
processing  equipment.  The  individual  devices  in  the  IDECS  System  are 
designed  primarily  for  versatility  of  operation  and  interconnection  to  allow 
evaluation  of  different  image  processing  techniques.  The  basic  capabilities 
of  the  IDECS  System  have  been  added  under  various  research  contracts 
with  the  principal  contributions  as  follows: 

1.  Real-time  video  processing  devices,  color  matrixing,  and 
initial  research  were  supported  primarily  under  NASA  contract 
#17-004-003. 

2.  A  six  image  synchronous  flying  spot  scanner  system  and 
general  purpose  real-time  decision  logic  with  color  coding 
capabilities  (the  spectrance  selector)  are  being  supported 
under  USAETL  contract  #  DAAK02-67-C-0435 . 

3.  Improved  operator  control,  digital  interface  and  the  capabilities 
for  automatic  control  card  processing  are  being  supported  under 
project  Themis . 

The  final  configuration  of  the  IDECS  System  as  it  is  now  envisaged  is  shown 
in  Figure  1 .  The  subsystems  which  are  being  supported  by  Project  Themis 
are  indicated  by  asterisks  and  are  described  individually  in  the  text.  The 
function  of  each  subsystem  is  the  same  as  that  presented  in  the  semiannual 
report;  however,  the  status  of  each  is  different. 

I.  ANALOG -TO -DIGITAL  CONVERSION  SUBSYSTEM 

A.  Function.  The  function  of  this  subsystem  is  to  convert  analog 
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signals  from  the  scanning  densitometer  of  any  of  the  flying -spot  scanners 
into  digital  information.  Each  analog  value  is  converted  to  a  binary  number 
variable  between  six  and  twelve  binary  digits  and  then  stored  in  one  or  two 
6-bit  words  on  magnetic  tape.  This  subsystem  is  also  to  have  sufficient 
flexibility  so  that  several  channels  can  be  multiplexed  together.  This  would 
allow  the  repetitive  sampling  of  each  flying-spot  scanner  in  a  sequence. 

B.  S ta tus .  The  analog  to  digital  converter  and  the  magnetic  tape 
unit  have  been  received;  however,  there  was  considerably  delay  because 
the  magnetic  tape  unit  had  to  be  returned  to  the  factory  for  repairs.  The 
design  for  the  high  speed  mode  has  been  completed  and  the  printed  circuit 
boards  are  being  constructed.  The  design  for  the  increment  mode  is  nearing 
completion.  The  total  interface  circuitry  for  the  analog  to  digital  conversion 
should  be  finished  shortly  after  the  cabinet  enclosure  arrives  which  will  be 
about  thirty  days  after  the  money  for  the  next  fiscal  year  is  received.  A 
block  diagram  of  this  subsystem  is  shown  in  Figure  2 . 

II.  DIGITAL-TO-ANALOG  SUBSYSTEM 

A.  Function .  The  function  of  this  subsystem  is  to  convert  digital 
information  from  a  digital  magnetic  tape,  a  small  computer,  or  other  digital 
information  sources  into  an  analog  voltage.  Initially  this  subsystem  will 
be  able  to  convert  two  binary  numbers  of  ten  bits  each  into  two  analog  vol¬ 
tages  suitable  for  input  into  a  CRT  display. 

B.  Status .  The  two  required  digital-to-analog  modules  have  been 
ordered  and  received.  Some  necessary  control  circuits  have  been  designed 
and  implemented;  however,  there  still  is  a  considerable  amount  of  design 
remaining  to  be  completed.  The  block  diagram  for  this  subsystem  is  shown 
in  Figure  3 . 

III.  DISPLAY  INTERFACE  SUBSYSTEM 

A.  Function .  The  function  of  this  subsystem  is  to  enter  data  at  a 
slow  rate  into  the  disk  memory  and  then  display  this  information  at  video 
rates  onto  the  color  display  unit  or  any  black-and-white  monitor.  Also 
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the  subsystem  will  eventually  have  the  capability  ot  being  usod  (or  data 
storage  in  a  small  computer. 

B.  Status .  The  disk  momory  has  boen  received;  however,  some  of 
the  prevloui  doslgn  objectives  had  to  be  altered  in  ord'i  to  obtain  suitable 
reliability  In  the  vidoo  storage  unit.  A  tentative  design  (or  tho  signal  con¬ 
ditioning  circuits  betweon  the  dlsK  mommy  and  the  vidoo  disolay  has  boen 
completed;  however,  there  is  still  considerable  design  work  remaining  on 
the  control  circuitry.  This  subsystem  Is  shown  in  Figure  4. 

IV.  SMALL  CQMfVIER  INIfiREACR  gUMXSIEM 

A.  Function .  The  purpose  of  tnls  subsystem  is  to  provide  a  computer 
Interface  so  that  data  may  be  processed  and  various  control  functions  per¬ 
formed  in  "real  time"  or  noar  "real  time." 

B.  Status .  The  effort  on  this  subsystem  will  be  started  when  tho 
necessary  funds  for  a  small  computer  are  obtained. 

V.  SCANNING  PENSlIQMSlSfl 

A.  Function .  To  process  an  Image  on  a  digital  computer.  It  Is 
required  that  the  Image  Information  be  changed  to  an  analog  signal.  This 
can  be  done  by  using  a  scanning  densitometer  which  consists  of  a  rotating 
and  translating  glass  drum  with  a  photomultiplier  tube  Inside.  The  Image 
is  placed  around  the  glass  drum.  A  light  source  placed  above  the  drum 
shines  through  the  image  Into  the  photomultiplier  tube.  The  output  of  the 
photomultiplier  is  an  analog  signal  proportional  to  the  grey  scale  of  the 
linage.  The  translation  of  the  Image  and  the  spot  size  are  set  so  that  20 
line  pairs/mm  are  produced . 

B.  Status .  The  basic  unit  has  been  designed  and  built;  however, 
thore  is  currently  some  difficulty  with  a  photo-diode  which  will  be  replaced 
as  soon  as  funds  become  available. 

VI.  MASJELYIPEQ.CPKIRQL 

A.  Function .  To  evaluate  the  functions  and  subsystems  individually 
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and  combinationally ,  it  is  nocossory  that  the  IDCCS  system  have  the  capabi¬ 
lity  to  connoct  and  disconnect  any  device.  This  is  the  (unction  of  the  vldoo 
control.  By  switching  the  vidoo  control,  various  combinations  of  input  and 
processing  devices  are  presented  at  the  cutput  devices.  The  ID£CS  block 
diagram  (Figure  1)  shows  the  input,  processing,  and  output  devices  thatare 
connoctod  to  tho  video  control. 

B.  Status.  Tho  vldoo  control  has  been  designed,  constructed  and 
connoctod  into  tho  system. 

VII.  POSITION  AND  SIZE  DEVICES 

A.  Function.  In  order  to  combine  and  process  more  than  ono  imago, 
it  is  necessary  that  the  lmagos  be  geometrically  congruenced.  This  can  be 
accomplished  by  controlling  the  raster  on  the  flying  spot  scanning  units. 

Tho  position  and  size  devices  control  this  raster. 

Tho  positioning  device  will  translate  and  rotate  tho  roster,  causing 
tho  image  on  the  display  to  movo  horizontally,  vertically,  and  angularly. 

Tho  size  dovlco  will  expand  or  contract  tho  raster;  therefore,  the  size  of 
the  imago  on  tho  display  is  controlled  horizontally  and  vertically.  The 
position  and  size  dovlcos  are  shown  in  Figure  5. 

This  'Jovlco  can  accomodate  the  following  corrections  for  a  8  x  10  cm. 
image  field: 

1.  Vertical  position:  +  3  cm. 

2.  Horizontal  position:  +  5  cm. 

3.  Vertical  size;  3  :  1 

4.  Horizontal  size:  3  :  1 

5.  Rotation:  +15* 

B.  Status.  Both  the  position  and  size  devices  have  been  imp’e- 
mented  for  ell  six  flying  spot  scanning  units. 

VIII.  SYNCHRONIZATION  DEVICE 

A.  Function  .  For  the  information  on  the  input  and  output  devices  to 
correspond,  it  is  essential  that  these  devices  scan  in  unison.  The  synchro¬ 
nization  device  accomplishes  this  by  sending  pulses  to  the  horizontal  input 
of  the  sweep  circuits  at  a  rate  of  15,750  Hz  and  to  the  vertical  input  of  the 


sweep  circuits  at  a  rate  of  60  Hz.  The  synchronization  device  Is  shown  In 
Figure  S.  This  device  Is  crystal  controlled  to  a  frequency  accu  acy  of 
.003  %. 

B.  Status.  A  synchronization  generator  Is  In  operation  to  match  the 
flying  spot  scanning  units  to  the  displays. 

DC.  fiQLfit  miLfiL 

A.  Function.  The  color  display  Is  i  sed  to  convert  video  signals 
proportional  to  gray  scale  on  photographs  or  -adar  Images  Into  a  color  which 
varies  In  Intensity  and  hue. 

B.  Status.  The  color  display  Is  operational. 

X.  TEXTURE  GENERATOR 

A.  Function  .  The  texture  generator  produces  an  electronic  equi¬ 
valent  of  alp-a-tone  patterns.  This  allows  multiple  categories  to  be  coded 
In  gray  levels  and  patterns  rather  than  color. 

B.  Status.  Although,  no  funding  has  been  made  available  to  build 
the  texture  generator,  Initial  design  and  feasibility  studies  have  been 
completed  with  in-house  funds. 

XT.  IBCmRS-EILTlSB- 

A.  Function  .  The  texture  filter  Is  used  to  discriminate  Image  filters 
on  the  basis  of  Image  textures  rather  than  Image  gray  levels. 

B.  Status.  Feasibility  and  design  for  the  texture  filters  are  complete 

from  in-house  funds.  A  prototype  model  Is  being  constructed  and  tested  under 
* 

a  NASA  contract. 
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Figure  2.  ANALOG  TO  DIGITAL  CONVERSION  SUBSYSTEM 
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Figure  5.  POSITION,  SIZE,  AND  SYNCHRONIZATION  DEVICES 
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for  the  Discrimination  of  Terrain  Characteristics,  "  R.  D. 
Ellermeier,  D.  S.  Simonett,  L.  F.  Dellwig.  Presented  at 
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lished  in  the  1967  IEEE  International  Convention  Record,  Part 
2,  p.  128-135.  Supported  by  NASA  Contract  NSR  17-004-003  and 
NASA  Grant  NsG-298. 

CRES  Technical  Report  61-20  "Evaluation  of  Multiple  Polarized  Radar 
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NSR  17-004-003  and  NAS  9-7175. 

CRES  Technical  Report  61-21  "The  Potential  of  Radar  as  a  Remote  Sensor 
in  Agriculture:  1.  A  Study  with  K-Band  Imagery  in  Western 
Kansas,"  D.  S.  Simonett,  J.  E.  Eagleman,  A.  B,  Erhart,  D. 

C.  Rhodes  and  D.  E.  Schwarz,  March  1967.  Supported  by 
NASA  Contract  NSR  17-004-003  and  NASA  Grant  NsG-298. 

CRES  Technical  Report  61-22  "Field  Studies  on  Vegetation  at  Horsefly 
Mountain,  Oregon  and  Its  Relation  to  Radar  Imagery,"  S.  A. 
Morain,  January  1967.  Supported  by  NASA  Contract  NSR 
17-004-003  and  NASA  Grant  NsG-298. 

CRES  Technical  Report  61-23  "K-Band  Radar  in  Vegetation  Mapping, " 

S.  A.  Morain  and  D.  S.  Simonett.  Presented  at  the  Annual 
Convention  of  the  American  Society  of  Photogramme  try,  Wash¬ 
ington,  D.  C.,  March  1967.  Published  in  Photogrammetrlc 
Engineering,  vol.  33,  no.  7,  July  1967,  p.  730-740. 

Supported  by  NASA  Contract  NSR  17-004-003  and  NASA  Grant 
NsG-298. 

CRES  Technical  Report  61-24  "Geologic  Evaluation  by  Radar  of  NASA 

Sedimentary  Test  Site,"  H.  C.  MacDonald,  P.  A.  Brennan, 
and  L.  F.  Dellwig.  Published  in  IEEE  Transactions  on  Geo¬ 
science  Electronics,  vol.  GE-5,  no.  3,  December  1967, 
p.  72-78.  Supported  by  NASA  Contract  NSR  17-004-003  and 
NASA  Grant  NsG-298. 
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NASA  Contract  NSR  17-004-003. 
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CRES  Technical  Report  61-27  "Application  of  Color- Combined  Multiple- 
Polarization  Radar  Images  to  Geoscience  Problems,"  D.  S. 
Simonett,  December  1966 .  Presented  at  the  Colloquium  on 
Classification  Procedures:  Computer  Applications  in  the  Earth 
Sciences,  Kansas  Geological  Survey,  Lawrence,  Kansas, 
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lished  as  Computer  Contribution  7,  State  Geological  Survey, 
University  of  Kansas,  Lawrence,  Kansas,  December  1966, 
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Huntsville,  Alabama,  June  1967.  Published  in  the  Transactions, 
vol •  II,  p.  MM- 23-1  to  MM-23-20 ,  Supported  by  NASA  Contract 
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Supported  by  NASA  Contract  NSR  17-004-003,  NASA  Grant 
NsG-298,  and  NASA  Contract  NAS  9-7175. 
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Paper  no  67-767,  21  pp.  Supported  by  NASA  Contract  NSR 
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P.  Youngberg,  August  1967.  Supported  by  NASA  Contract  NAS  9-7175. 
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CRES  Technical  Memorandum  61-3  "A  Method  of  Analyzing  NRL  Flight 

Radar  Return,"  J.  W.  Rouse,  Jr.,  September  1964.  Supported 
by  NASA  Contract  NSR  17-004-003. 
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October  1964.  Supported  by  NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-5  "  Initial  Design  Considerations  for  the 
Multi-Color  Display,"  G.  W.  Dalke,  J.  E.  Rathke,  J.  W. 

Rouse,  Jr. ,  December  1964.  Supported  by  NASA  Contract  NSR 
17-004-003. 

CRES  Technical  Memorandum  61-6  "Establishment  of  the  Criteria  for 
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Beam  Width  and  Pulse  Length  on  the  Ground,"  J.  Butler, 

October  1964.  Supported  by  NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-8  "Synthetic  Aperture  Power  Requirements," 
R.  K.  Moore,  J.  Butler,  November  1964.  Supported  by  NASA 
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J.  Butler,  R.  K.  Moore,  March  1965.  Supported  by  NASA 
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NSR  17-004-003. 

CRES  Technical  Memorandum  61-13  "Synthetic  Aperture  Antenna,"  J.  W. 
Rouse,  Jr. ,  May  1965.  Supported  by  NASA  Contract  NSR 
17-004-003. 

CRES  Technical  Memorandum  61-14  "Calculation  of  Tristimulus  Coefficients 
for  the  CRES  Multicolor  Display,"  G.  W.  Dalke,  May  1965. 
Supported  by  NASA  Contract  NSR  17-004-003. 
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CRES  Technical  Memorandum  61-15  11  Modifications  Made  on  Photovolt  520-A 
Densitometer,"  M.  Buchanan,  May  1965.  Supported  by  NASA 
Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-16  "Unfocused  Synthetic  Aperture  Antenna 
Radar  Power  Requirements , "  J.  W.  Rouse,  Jr.,  June  1965. 
Supported  by  NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-17  "Theory  of  Synthetic  Aperture," 

M.  Buchanan,  June  1965.  Supported  by  NASA  Contract 
NSR  17-004-003. 

CRES  Technical  Memorandum  61-18  "Remote  Sensor  Studies  of  the  Pisgah 
Crater  Area,  California:  A  Preliminary  Report, "  L.  F.  Dellwig, 

M.  E.  Bickford,  N.  Kirk,  R.  Walters,  December  1965 . 

Supported  by  NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-19  "Technique  for  Producing  a  Pseudo 
Three-Dimensional  Effect  with  Monoscopic  Radar  Imagery," 

L.  F.  Dellwig,  H.  C.  MacDonald,  J.  N.  Kirk,  December  1965. 
Supported  by  NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-21  "Reduction  of  Like-Images  to  Space 
Correlated  Data  Lists,"  G.  W.  Dalke,  M.  O.  Vichellich, 

May  1966.  Supported  by  NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-22  "The  [3x3]  Matrix  Unit, "  F.  Welter, 
May  1966.  Supported  by  NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-23  "Pattern  Recognition  Problem,"  R.  M. 

Haralicfr,  June  1966.  Supported  by  NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-24  "Photographic  Color  Combination 
Procedures  and  Equipment  Available  at  CRES,"  D.  Egbert, 

June  1966.  Supported  by  NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-25  "Radar  Bibliography  for  Geoscientists," 

R.  L.  Walters,  June  1966.  Supported  by  NASA  Contract  NSR 
17-004-003. 

CRES  Technical  Memorandum  61-26  "Radar  Linears  not  Related  to  Geology," 

H.  D.  MacDonald,  E.  Gillerman,  L.  F.  Dellwig,  J.  N.  Kirk, 

July  1966.  Supported  by  NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-27  "Radar  Fundamentals, "  J.  W.  Rouse, 

Jr.,  W.  P.  Waite  (being  completed).  Supported  by  NASA 
Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-28  "Image  Reproduction  with  Color  Tele¬ 
vision  Receiver,"  R.  L.  Knuckey,  July  1966.  Supported  by 
NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-30  (Deleted) 
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CRES  Technical  Memorandum  61-31  "Automatic  Acuteness  Measurements 
for  Multiple  Images,"  G.  W.  Dalke,  August  1966.  Supported 
by  NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-32  "Automatic  Measurement  and  Reduc¬ 
tion  of  Drainage  Basin  Data,"  G.  W.  Dalke,  August  1966. 

Supported  by  NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-33  "Automatic  Measurements  of  Texture 
and  Detection  of  Lines  in  a  Flying- Spot  Scanner  System," 

G.  W.  Dalke,  August  1966.  Supported  by  NASA  Contract 
NSR  17-004-003. 

CRES  Technical  Memorandum  61-34  (Deleted) 

CRES  Technical  Memorandum  61-35  "Pattern  Recognition  and  Similarity," 

R.  M.  Haralick,  August  1966.  Supported  by  NASA  Contract 
NSR  17-004-003. 

CRES  Technical  Memorandum  61-36  "Design  System  for  Flying  -  Spot 
Scanner,"  J.  E.  Rathke,  August  1966.  Supported  by  NASA 
Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-37  "The  CRES  Data  Handling  System," 

J.  E.  Rathke,  August  1966.  Supported  by  NASA  Contract 
NSR  17-004-003. 

CRES  Technical  Memorandum  61-38  "Photomultiplier,  Compensated 

Video  Amplifier  and  Video  Amplifier  Circuits,"  J.  £.  Rathke, 
September  1966.  Supported  by  NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-39  "The  3-Channel  Flying-Spot  Scanner, 

J.  E.  Rathke,  August  1966.  Supported  by  NASA  Contract  NSR 
17-004-003. 

CRES  Technical  Memorandum  61-40  "The  Tri-Color  Display  Unit,"  J.  E. 

Rathke,  September  1966.  Supported  by  NASA  Contract  NSR 
17-004-003. 

CRES  Technical  Memorandum  61-41  (Deleted) 

CRES  Technical  Memorandum  61-42  "Multi-Color  Presentation  of  Images 
with  Edge  Enhancement,"  R.  L.  Knuckey,  October  1966. 

Supported  by  NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-43  "Synthetic- Aperture  Signal  Film  Quantity 

Calculations,"  J.  W.  Rouse,  Jr.,  October  1966.  Supported  by  NASA 
Contract  NSR  17-004-003. 
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ORES  Technical  Memorandum  61-44  "An  Analog  to  Digital  to  Paper  Tape 
Converter,"  J.  Betlack,  April  1966.  Supported  by  NASA 
Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-45  "Terrain  Discrimination  by  Radar  Image 
Polarization  Comparison,"  R.  K.  Moore,  L.  F.  Dellwig,  October 
1966.  Published  in  Proceedings  IEEE,  vol.  54,  no.  9,  September 
1966,  p.  1213-1214.  Supported  by  NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-46  "Analog  to  Digital  Converter  for  Low 
Frequency  Sampling,"  R.  Field,  April  1967.  Supported  by 
NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-47  "Evaluation  of  Geoscience  Information 
as  a  Function  of  Radar  Resolution,"  R.  McCoy  (Classified- 
internal  distribution  only) .  Supported  by  NASA  Contract  NSR 
17-004-003. 


CRES  Technical  Memorandum  61-48  "Electronic  Processing  for  Synthetic 
Aperture  Array  Radar,"  R.  K.  Moore,  J.  W.  Rouse,  Jr., 
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p.  2  33-234.  Supported  by  NASA  Contract  NSR  17-004-003. 

CRES  Technical  Memorandum  61-49  (Deleted) 

CRES  Technical  Memorandum  61-50  "Radar  Imagery,  a  New  Tool  for  the 
Geologist,"  J.  W.  Kirk,  R.  L.  Walters.  Published  in  The 
Compass  of  Sigma  Gamma  Epsilon,  vol.  43,  no.  2,  January 
1366,  p.  85-93.  Supported  by  NASA  Contract  NSR  17-004-003. 
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Supported  by  NASA  Contracts 
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P.  Youngberg,  July  1967.  Supported  by  NASA  Contracts  NSR 
17-004-003  and  NAS  9-7175. 

CRES  Technical  Memorandum  61-56  "IBM  Paper  Tape  Punch  Program  for 

the  LGF-30  Computer,"  P.  Youngberg,  August  1967.  Supported 
by  NASA  Contracts  NSR  17-004-003  and  NAS  9-7175. 
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Printer  for  the  LGP-30  Computer,"  P.  Youngberg,  August  1967. 
Supported  by  NASA  Contracts  NSR  17-004-003  and  NAS  9-7175. 
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NASA/MSC  Mission  77,  Site  76,  Garden  City,  Kansas,"  J.  Ratzlaff, 
October  1968.  Supported  by  U.S.G.S.  Contract  No.  14-08-0001-10848. 
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CRES  Technical  Memorandum  118-1  "Pattern  Recognition  of  Multi -Spectral 

Photographs  Using  a  Conditional  Probability  Method, "  R.  Haralick, 
September  1967.  Supported  by  NASA  Contract  NAS  9-7175. 

CRES  Technical  Memorandum  118-2  "Line  Driving  Amplifier,"  J.  Bryant, 

November  1967.  Supported  by  NASA  Contract  NAS  9-7175. 

CRES  Technical  Memorandum  118-3  "Synthetic  Aperture  Radar  Acoustic 

Simulator,"  L.  Fromme,  December  1967.  Supported  by  NASA 
Contract  NAS  9-7175. 

CRES  Technical  Memorandum  118-4  "Preliminary  Mission  Analysis  Report; 

NASA/MSC  Mission  No.  47,  Site  93,  Pt.  Barrow,  Alaska," 

J.  Rouse,  Jr.,  October,  1967.  Supported  by  NASA  Contract  NAS  9- 
7175. 

CRES  Technical  Memorandum  118-5  "Preliminary  Report  on  Radar  Lineaments 
in  the  Bosto  Mountains  of  Arkansas,  "  J.  Kirk  and  R.  Walters. 
Published  in  The  Compass  of  Slama  Gamma  Epsilon,  vol.  45, 
no.  2,  January  1968,  p.  122-127.  Supported  by  1  ASA  Contract 
NSR  17-004-003  and  NASA  Grant  NsG-298. 

CRES  Technical  Memorandum  118-6  "On  the  Feasibility  of  Imaging  Radar  for 
Small  Spacecraft,  "  J.  Rouse,  Jr.,  (being  completed) .  Supported 
by  NASA  Contract  NAS  9-7175. 

CRES  Technical  Memorandum  118-7  "Use  of  the  Radar  Scatterometer  in  Aircraft 
and  Spacecraft  Resources  Programs,"  R.  K.  Moore,  October  1968. 
Supported  by  NASA  Contract  NAS  9-7175. 

CRES  Technical  Memorandum  118-8  "National  Aeronautics  and  Space  Administra¬ 
tion  Earth  Resources  Survey  Program  90-Day  Mission  Analysis  Report, 
NASA/MSC  Mission  70,  Site  166,  Iceland,"  W.  P.  Waite  and  R.  L. 
Walters,  August  12,  1968.  Supported  by  NASA  Contract  No.  NAS  9-7175 
and  Naval  Oceanographic  Contract  N62306-67C-0044. 

CRES  Technical  Memorandum  118-9  "National  Aeronautics  and  Space  Administra¬ 
tion  Earth  Resources  Suivey  Program  90-Day  Mission  Analysis  Report, 
NASA/MSC  Mission  74,  Site  76,  Garden  City,  Kansas,"  D.  H.  Berger, 
17  September  1968.  Supported  by  NASA  Contract  NAS  9-7175  and 
U.S.G.S.  Contract  14-08-0001-10848. 

CRES  Technical  Memorandum  118-10  "National  Aeronautics  and  Space  Administra¬ 
tion  Earth  Resrouces  Survey  Program  90-Day  Mission  Analysis  Report, 
NASA/MSC  Mission  74,  Site  85,  Lawrence,  Kansas,"  R.  L.  Walters, 
September  1968.  Supported  by  NASA  Contract  NAS  9-7175  and  U.S.G.S. 
Contract  No.  14-08-0001-10848. 
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CRES  Technical  Memorandum  118-11  "National  Aeronautics  and  Space  Administra¬ 
tion  Earth  Resources  Survey  Program  90-Day  Mission  Analysis  Report, 
NASA/MSC  Mission  74,  Site  85,  Lawrence,  Kansas,"  R.  L.  Walters, 
September  1968.  Supported  by  NASA  Contract  No.  NAS  9-7175  and 
U.S.G.S.  Contract  No.  14-08-0001-10848. 

CRES  Technical  Memorandum  118-12  "National  Aeronautics  and  Space  Administra¬ 
tion  Earth  Resrouces  Survey  Program  90-Day  Mission  Analysis  Report, 
NASA/MSC  Mission  54,  Site  85,  Lawrence,  Kansas,"  W.  G.  crooner 
and  R.  L.  Walters,  (In  preparation).  Supported  by  NASA  Contract 
No.  NAS  9-7175  and  U.S.G.S.  Contract  No.  14-08-0001-10848. 
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CRES  Technical  Memorandum  112-1  "Acoustic  Scattering  from  Surfaces  with 
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